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PREFACE 


This book contains the substance of a course of lectures 
delivered to engineers and physicists at the National 
Standards Laboratory, Sydney, in 1944, when the author 
was seconded to the Radiophysics Laboratory of the 
Council for Scientific and Industrial Research. A mimeo- 
graphed edition was produced by the C.S.LR. in 1945 
for circulation in Australia. 

There is an increasing demand for a knowledge of 
operational methods by engineers and physicists who 
encounter transient problems, but whose mathematical 
equipment is restricted to that of an ordinary engineering 
course. Of the many books available, the older ones use 
the rather suspect and certainly obscure “ operational 
methods while the newer ones require a knowledge of 
the theory of functions of a complex variable, and some 
comparatively advanced mathematics. Here I have at- 
tempted to go as far as possible using no mathematics 
beyond ordinary calculus, but at the same time catering for 
the reader who wishes to acquire the manipulative skill 
necessary to solve his own problems : to this end, the book 
contains as little theory as possible ; it is, in fact, largely a 
collection of worked examples illustrating the methods of 
solution of the various types of problem commonly arising 
in circuit theory ; a handful of problems, with answers, 
has been added to enable the reader to test his skill. 

For problems on ordinary linear differential equations, 
and thus on electric circuits with lumped constants, the 
use of the theory of functions of a complex variable is 
quite unnecessary ; but it is needed for a complete study 
of partial differential equations, and in the higher parts 
of the subject. The engineer who wishes to go deeply 
into these matters must study this theory, and since quite 
sophisticated points in it do arise, he must study it fully 
and carefully, and so before he starts he is entitled to 
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know what types of problem he will be able to solve, and 
how useful the solutions will be. I have therefore ^ m 
Chapter IV discussed the equations of the transmission 
line from the point of view of the earlier chapters, using 
calculus only ; this treatment covers all the common 
problems, and will enable the reader to follow published 
work on the subject. It is to be understood that solutions 
obtained in this way are correct, and that the algebra is 
the same as it would be if the complex variable were used : 
the latter merely changes the point of view, and allows n 
complete mathematical discussion of the validity of the 
solutions. 

The reasons for the notation used here for the Laplace 
transform ^oo , v 

\e-t>*y{t)dt . . (i) 

of 3'(<) may be given briefly. Any notation should satisfy 
three conditions : (i) it must be short, and easy to write J 
(ii) it must be applicable to any symbols ; (iii) it should 
not be used in any other connection. The simplest 
notation of all' would be to write a capital letter Y for the 
Laplace transform of the corresponding small letter y ; 
but this may be most inconvenient since it is not applicable 
to many symbols, in particular to capital and many Oreek 
letters. A “ star y*, is not much used in other con- 
nections but is rather clumsy to write. The "‘bar 
notation, jJ, used here, was originally used by Ldvy and 
reintroduced by Carslaw ; it is easy to write, and can 
be applied to any symbols, but it has the disadvantage of 
being employed in other senses, these, however, rarely 
occur in those fields where the Laplace transformation 
is most useful, and ease of writing is a most importanl: 
desideratum. If preferred it can be replaced by y used 
by Millman. 

A much more difficult question is that of the relative 
merits of the form (i) and of the “^-multiplied ” form 

fOO 

p\ e~^^y(t)dt , . . . fsa j 
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Both notations have their own minor advantages, but the 
essential point is that form (2) gives an exact correspondence 
with the large body of work done by the older operational 
methods, while (i) corresponds with the newer books, 
with most of the theoretical work, and with corresponding 
work on the Fourier transform. 

Finally, it is a pleasure to acknowledge the assistance 
given to me by many officers of the C.S.I.R., in particular 
that of the Secretary in making available the drawings 
for the text-figures. 
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CHAPTER I 


FUNDAMENTAL THEORY 
I* The Laplace transform 

In this book we shall primarily be concerned with 
problems on the behaviour of 'dynamical or electrical 
systems which are started in a given way at some instant, 
which we take to be the zero of time, i = o. Mathemati- 
cally, these are called ‘‘ initial value problems ”, and may be 
expressed as the solution of a system of differential equa- 
tions with given initial values at the instant t = o. Negative 
values of the time do not arise in either the specification of 
the problems or in their solution. The Laplace trans- 
formation is a mathematical device which is useful for 
solving such problems : if a function is specified for all 
positive values of the time, we can (in the simple cases 
needed here) write down a related quantity called its 
Laplace transform ; conversely, if the Laplace transform 
is known, the value of the function can be found from it. 
The importance of the method arises from the fact that, 
for the types of problem we are interested in, it is easy to 
write down the Laplace transform of the solution, and 
from this the solution itself is then easily found. 

We suppose y{t) to be a known function of t for values 
of t> 0. Then the Laplace transform y{p) of y{t) is 
defined as 

y(P) = f . • (i-i) 

where is a number ♦ sufficiently large to make the integral 
(i.i) convergent. For example, if y{t) = i any positive 

* Of course p may be a complex number whose real part is 
sufficiently large to make (i.i) convergent, but, in the early parts 
of the theory, it is a little more definite to think of p as a real 
positive number, and nothing is lost by this restriction. The 
Laplace transform is one of a number of integral transforms 
formed as in (i.i) by multiplying a function of £ by a chosen 

I 
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p is permissible, but if y{t) = we must have p > 
Apart from this p is unrestricted, and y{p) is a function 
of p just as y{t) is of L We shall occasionally write them 
in full in this way to emphasize this dependence, but 
usually shall write y for the Laplace transform of 3^, it 
being understood that these. are functions of p and t 
respectively. 

We need first a collection of the Laplace transforms 
of a few common functions. For example, if jy = 
where a may be real or complex, 

TOO fOO I 

y = 1 . e^^dt = I = , 

Jo Jo P ^ 


where p must be greater than the real part of a. 
Again, if a is real and 

y = sm at — 
using (1.2) we have 

I / I 


zi\p —ia p p^ + d 


( 1 * 3 ) 


Proceeding in this way, or quoting the results from 
known definite integrals, we can build up the Table I 
of Laplace transforms which contains all the transforms 
needed in the first three chapters. This Table, and 
Theorems I to III below, provide the whole of the theory 
needed in Chapters I and II. 

In (1.6), as remarked above, p has to be greater than 
the real part of a, while in (1.7) to {1.12) a is real, and in 
(1.9) and (i.io) p must be greater than* \a\. All the 
results of Table I are elementary integrals except (i.il) 
and (i.12), which have been included because they arise 
in cases of resonance. The simplest way of deriving 


function of t and py and integrating with respect to t* The 
Fourier transform (26.1), is another example. Each integral 
transform has characteristic properties, but those of the Laplace 
transfoim which are proved below make it the most suitable for 
the solution of initial value problems. 

* If a is real, |fl| = a, if a > o ; ' |a| — — iz, if ^ < o. 
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results of this type is as follows : written out in full (i.8) 
states that 

poo jy 

1 cos at dt =: 

Jo p^-^a- 


Table I 


y{t) 

yiP) 


I 

1 

(1-4) 

(n - i) 1 

I ... 

— , n a positive integer 

pn 

(1.5) 


I 

p — a 

. (1.6) 

sin at 

a 

p^+ ' 

(i.?) 

cos at 1 

P 

p'^ a^ 

(i.8> 

sinh at 

a 

p'^ - • 

(i.9) 

t 

cosh at 

P . 

(i.io) 


p" — 

i ^ 

— sin at 

P 

(i.ii) 

2a 

(p- -h a'^y- 

-^-3 (sin at — at cos at) 

1 

(p^ + ay ■ 

(i.ia) 


If we differentiate * both sides of fhis with respect to Uy 
as if a were a variable, we get 


e-’PH sin at dt ■ 


zap 


ip-^ + a'Y 


* This technique of differentiating under the integral sign 
with respect to a parameter is a powerful weapon of pure mathe- 
matics ; naturally its use needs justification, and this can usually 
be supplied when, as here, it is employed in deducing new Laplace 
transforms from known ones. 
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which is eqmvalent to (i.ii). Treating (1.7) in the same 
way gives (1.12). 

The collection of transforms in Table I can be greatly* 
extended by the use of the following Theorem. 

Theorem I. If y{p) is the Laplace transform of 
and a is any number y real or complex, then y(p + a) is the 
Laplace transform of e’'°‘^y{t). 

This follows immediately, since the Laplace transform 
of e^^^y{t) is 

I . e--^*y(t)dt = | e-^^^-^°'^^y(t)dt — y{p + a). 

As an example of the use of this Theorem, we see from 
(1.5) that the transform of 

^-at ffi-l j 

(« — i)! (P + «)"' 

Similarly from Theorem I and (1.8) it follows that the 
transform of 

. . ' . p b 

cos at IS . , r^<. ■ — v 

{p + by + a^ 

and from Theoreln I and (1.7) that the transform of 


6“^* sin at is 


a 

(P + 


2. To find the function y{t) which has a given Laplace 
transform. 

It will be seen in §4 that the Laplace transformation 
method of solving an ordinary linear differential equation 
with constant coefficients and given initial conditions 
consists of writing down an algebraic equation for the 
transform y{p) of the solution y[t). This transform will 
usually be found to have the form 


fM 

iipi 


(2.1) 


where f(p) and g(p) are polynomials in p which have no 
common factor, the degree oif{p) being lower than that of 
g(.p)- 
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Most of the labour * of the solution consists of finding 
the y{t) corresponding to a y{p) given in this form. This 
is done by looking up y{p) in the Table I of transforms, 
possibly combined with Theorem I. It is shown at the 
end of this secton that the yip) found in this way is the 
only function which has y{p) for transform. 

For example, if 




I 


it follows from (i.6) that 3; = 


Similarly, from (1.7) and (1.8), if 

_ + 6 

y = p y = ZCOS2t + 2 sm 2 t. 

If the denominator of is a general quadratic in p, 
we complete the square, and use Theorem I and the Table. 
For example, if 

P + 7 ^ P + 7 ^ P + i , 6 

^ P^ + 2>p + 5 (P+i)^ + 4 (^+i)^ + 4'^(^+i)" + 4’ 

it follows from Theorem I and the Table that 
y = cos zt + 36”^ sin zt. 


Notice that both numerator and denominator of y have to 
be expressed in terms of + a (in the case above, ^ + i) 
before Theorem I is used. 

If the denominator of y is of higher degree than the 
second, we have first to find its factors, f then we express 
y in partial fractions by the usual algebraic methods, 
and finally write down the function of t corresponding 
to each fraction using Table I. Alternatively, Theorem II 
below may be used when it is applicable. 


* The student is particularly advised to make himself skilful 
in the processes of this section. Labour can be saved by the 
use of a more extensive table of transforms than that given here, 
but transforms not included in the tables at present available 
frequently arise, and in any case a solid background of manipulative 
skill is extremely useful. 

t For the method of doing this in problems where they are 
not obvious see § 12. 
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^ P ^ 

Example I. 

Assume* ^ ^ =~+ 

p{p^ +4P + 8) p ^p^ + # + 8 

Then ^ +- 1 s A(p^ + + 8 ) + ^(Bjp + C) . (2.^) 

Equating coefficients of p^, py and the constant term, 
respectively, we have 

A + B = o, 

4A + C “ I) 

8 A = I. 

Solving these for A, B, and C, we find 

^ j>-4 _ _ I (j> + g)-6 , . 

y-%p 8(^*= + 4^ + 8) Zp 8[(^> + 2r + 4 r 

Then from Theorem I and the Table, we find 
y = ^ ^e”^*{cos 2 t — 2 sin 2^}. 

The process of equating coefficients used above is 
always available, but the results can often be obtained more 
quickly otherwise. For example, if we put ^ — o in the 
identity (2.2) we find A = 1/8. Then, putting A = 1/8 
in (2.2), this becomes BjJ) -|- C s -1(4 — p), and we get the 
result {2.3) obtamed before. 

Example 3. y = ^ 

. I _ A , B ,Cp + D 

^^’^^(i> + i)^(^^+4) i'+ I + 1)^^ />^ + 4 ‘ 

Then 1 ^ A(^+i)C^H4)+B(^'+4)+(C^+D)(/.+i)*. 
Equating coefficients we find 

A B = C = D = — 

* It is shown in text-books on algebra that if we assume a 
fraction with a constant numerator corresponding to each linear 
factor of the denominator, and a fraction with a linear numerator 
corresponding to each quadratic factor (or squared linear factor) 
of the denominator, and so on, we shall always have just the right 
number of equations to determine the unknown constants in 
these fractions. 
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Then, using the Table and Theorem I, we get 

y = — A cos sin %t. 

There is a simple general Theorem applicable to many 
important cases : 

Theorem ^ IL If y{p) where f{p) and g{p) 

are polynomials in p^ the degree off(p) being less than that of 
q{p), and if 

SiP) = (P— ^)(P - « 2 ) •••(/> - «„), ■ (2.4) 

where a^^ a^, . . a^ are constants^ which may he real or 
complex but must all be different, then 

-^a^^g'ictr) ■ ■ ■ ■ • • (“-S) 


_y I _ /(£r) 

^ p—c^T ' iar-aj) . . . {a,—ar.x){a*+aT-\) ■ ■ • (Or— “n) 


(2.6) 


_Y I r (f> - ^r)/(j>) i 
;^i>-ar‘L g(p) J,_„/ 


(a. 7 ) 


Also the function y(t) whose Laplace transform is y(p) is 


y(t)=y4r^‘^* 

fi 

=I 7 Z^ 


Kar)^’’^ 


“ (flr-ai) • ■ • . . . (a,-a„) 


(2.8) 

(2.9) 


glp) J 


e^r*. 

Jji-Or 


(2.10) 


(2.8) is the form taken by Heaviside^s Expansion 
Theorem in the present method.* Since (2.8) . . . (2.10) 
follow directly from (2.5) . . . (2.7), respectively, by 

* In Heaviside’s notation, and in the “^-multiplied ’’ form of 
the Laplace transformation, in which (in order to have an exact 
correspondence with the Heaviside notation) the Laplace trans- 
form is defined as p times the integral in (i.i), an additional 
term appears arising from this extra factor, p, and the simplicity 
of (2.8) is lost. 
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using (i.6), it appears that the essential content of TheoreiT^^ 
II is simply the formulae (2.5) . . . (27) for puttirtg 
f(p)/g(p) into partial fractions in the case in which g(p) 
no repeated factors. This is a well-known result,^ 
for completeness we give a proof below. 

All of (2.5) to (2.10) are useful ; (2.6) is easily remertx*^ 
bered by noticing that the coefficient of il(p — a^) is 
result of putting p ^ in f{p)lg{p)> omitting the facto** 
(P ~ ^r) of g(p) I (2*7) is simply another way of writing 
this statement. Similar remarks apply to (2.9) and (2.10)- 
The partial fraction formulae (2.5) and (2.6) are themselves 
very useful, cf. Example 4, below. With regard to (tz.Sy 
and (2.9), it is usually better to use (2.9) when g(p) has ^ 
few simple factors, as in the examples of this section- 
It sometimes happens that g(p) has complicated factors, 
which need not even be set out explicitly, as in the problems 
of § 1 6, and then (2.8) is the more useful. 

To derive (2.6) we put y(p) into partial fractions by tlae 
methods discussed earlier, by assuming 


f(p)^f A, 
Then Sd’) 


(z.ii) 


^ip) = '^K{p~‘h) • • • 

f-1 

Putting p = a^m this gives 

ai) . . . (ar-^S'r-OCar-^r+l) • • • («r— ««)» 

r = I, 2, . . Tit 

and substituting this value in (2.1 1) gives (2.6). To deduce 
(2.5), we notice that differentiating (2.4) with respect to 
p by the ordinary rule for differentiating a product gives 

• • • (P-«r-.l)0-«r+l) • • • (/►-««)• 

r-1 

Putting ^ in this, gives 

• • • {«r-«r-0(«r-«r+l) • • • (2- 1 a) 

Using (2.12) in (2.6) gives (2.5). 

• Cf, Gibson, Treatise on the Calculus ^ edn. 3 (1906), § 120. 
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Before proceeding further it is desirable to settle some 
questions of terminology. Firstly, suppose j?(p) is given 
in the form (2.1). The simplest case is that of (2.4) in which 
g{p) is a product of n different factors, 

g{p) = (p — a^)(p , (p ^ 


The numbers are the roots of the equation 

g{p) == o. They may also be spoken of as the zeros of 
g(p), since they are the values of p for which g{p) is zero. 
Finally, they are the values of p for which the function 
y{p) is infinite, and in this sense are described as the poles 
of y{p)* In the more general case in which g{p) has 
factors such as {p — a-^ is called a multiple (or more 

precisely an n-ple) root of g{p) = o, an ;z-ple zero of 
g{p), or a pole of order n of 

KP)^ 

Secondly, there is the ques- 
tion of the notation for a 
complex number. If 
= A? + fy 

is a complex number, we 
call X the real part of y 
the imaginary part of and 
X —iy the conjugate of z. 

The number may be repre- 
sented on the Argand diagram 
of Fig. I, and may also be 
specified by its modulus, |<s'| = the angle 

(f) defined * by any tzuo of 



Fig. I. 


^ = yy ^ = l • ■ ■ ( 2 - 13 ) 

We shall call this angle ^ the argument of z, written 
arg z ; other names for it are the amplitude of z, (amp 
or the phase of .s’. 

*'This definition still leaves (f) undetermined by 'an integral 
multiple of ztr ; this does not matter in the applications made 
here. If desired, </> can be restricted to the range tt > — tt, 

when it is called the Principal Value of arg z. 
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With the notation, (2.13), any complex number x 
may be written in the form 

» + • . (2.x 4 -) 

where cj) — arg (x + ty). 

The importance of this in the present connection is thsti: 
when y{t) is found from y(p) by the use of (a.9), it will 
frequently involve products or quotients of compl^^ 
numbers, and the most rapid way of reducing such 
expressions to a real form is by using (2.14) for each of 
the complex numbers (cf. Example 5, below). 

One pitfall of some importance must be mentioned t 
the definition of </► in (2* 13) is equivalent to * 

^ = arg(^ + iy) = tm-'^yjx), if > o\ - c- \ 

= TT + t2Ln'~\ylx), if X <oJ, 

thus, if we had written (j> == ta.n~\yjx) in (2*14) (and this 
is very frequently done), there might be an error of ^ 
in the argument of ty, leading to an error in sign 
the final result. This point is particularly important it% 
Tables of transforms, and in problems such as Example S 
below, where numerical values of some parameters are 
not specified and thus may be taken with either sign. 

We now give some examples on the use of Theorem II- 

ExAMPLE 3 . y=—, ^ 

p{p + ^){p+ ^){p + 3) 

By (3.6), 

I ■ I 

.2.3 .^’^(“-l). I.3.(P + I) 

■I- ^ 1 ' ^ 

(-2)(“-i) . I . (^+2) ^ (- 3 )(-~^)(“i)(/>+- 3 )^ 
and y follows from (1.6). Or, using (2,9), 

v= — ^ ^ g-2i 

^(-3)(-2)(-i)^ - 

* tan“^ X is defined as the angle between — J-nr and Jtt whose 
tangent is x. 
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Example a, y = - wy— o ; — o - . — v* 

^ p^ip^ + i){p^ + 4 ) 

Using (a.6) on j;, regarded as a function of we find 

y ~ 4p^ 2iP‘^ + i) i2(jf 2 -(- 4)- 

Thus V = — - sin i + — sin 2t 

4 3 24 

Problems on resistanceless circuits often lead to transforms 
which can be conveniently treated, as above, as functions 
of p^. 

Example 5. 

^ ^ + <2 

p 4 " ^ 

(p + p)ip + M — i^)ip + M + »»)• 


This is a generalization of Example i, and could be 
dealt with either by resolving into partial fractions with 
real denominators, as was done in Example i, or by re- 
solving into partial fractions with linear complex denomin- 
ators, that is, by using one of the forms of Theorem II. 
Thus applying (2.9) to (2. 16) gives * 


^ (f. - by + 


+ { (/_ + Conjugate} (2.17) 


Then writing 


4 i = arg[(a — /r) + in] 
= arg[(Z> — /a) + m]. 


* A pair of conjugate roots of g{p) = o give in (2.9) the sum of a 
pair of conjugate terms, of which only one need be written down. 
The method used below for reducing such a pair of terms to real 
form will be found very convenient. 
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and using (2.14) in the complex numbers in (2.17), we 
find 




ilJ^-by + n^^ 

I i( {a-tiy + n^ 

n\(6 — fjLy + 

(^^-by + n^^ 

, if (a-f<.y + n^ 
n\(b — /jt)2 -j- 


I 1_ ConjugateV 

I sin {nt + — ^2)- 


Theorem II covers the case * in which g{p) has no repeated 
factors ; for example, it is not applicable to Examples 2 
or 4 above. It is easy to generalize the theorem to the 
case of repeated factors. Theorem II may be stated in 
the form that, to each linear factor [p — af) of the denomin- 
ator of y{p) there corresponds a term 


, . . (2.18) 

. , S K) 

tn the solution. 

The generalization is that, to each squared factor 
{p —hy of the denominator of y(p) there corresponds a term 


[{P “ hff(p>)-\ 
L g{p) J 


te^^ + 

l )-6 


-d- ( {p-byf{p) -\- 
jp\ g{p) j_ 


3)“& 


in the solution. 


(2.19) 


* In (a.i) it was stated that/(^) and g{p) were to have no common 
factors : this is desirable but not essential. With complicated 
expressions (such as those of § 16) it would greatly increase the 
labour to separate out such a factor, and, in fact, if f{p) and g{p) 
both have the same factor {yp — a), the term in c®* corresponding 
to this in Theorena II has zero coefficient, and so the application 
of Theorem II gives the correct result. If g{p) contained a 
factor (p — a)®, and }{p) a factor {p — a), Theorem II could 
be applied after cancelling the factor — a), whereas a cursory 
inspection would suggest that it is not applicable and that (2.18) 
would have to be used (this, also, would give the correct result). 
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And to each m-ple factor (p — c)^ of the denominator 
of y{p) there corresponds a term 


r 

r“ 

1 


1 g{p) /j 


fm-s-l 


, 5 I (m — ^ — i) ! 


(2.20) 


in the solution. 

As an example, consider the problem of Example 2 
above, in which 


(p-\- +4) (i> + - 2i)[p + 2iy 

Then by (2.19) and (2.18) 

\.W+ 4 - 




leading to the same result as before. 

Finally, the question of '' uniqueness ** must be dis- 
cussed. To a given function y{t) there corresponds a 
single transform y{p) defined by (i.i). It is not clear, 
however, that to a given transform y{p) there corresponds 
only a single function y{t). Given y{p) we have sometimes 
been able to find a function y{t) which has j 7 (p) for transform 
by looking up y{p) in the Table of transforms ; but there 
is so far no reason why the function of t found this way 
should be the only function which has y{p) for transform. 
The same question arises in the integral calculus when 
we seek a function which has a given function for differential 
coefficient, and there it is found that there are infinitely 
'many such functions (differing by additive constants). 

In the present case the matter is settled definitely by 
Lerch's Theorem,^ which states that, if two continuous 


* For a proof see C. and J., Appendix I. The actual theorem 
ia a great deal wider than that stated above, and includes also 
the useful result that if two functions yi{t) and y^it) are continuous 
except at a number of ordinary discontinuities (i.e. jumps in value), 
and have the same Laplace transform, they must have the same 
points of discontinuity and can differ only at these points. 
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functions ham the same Laplace transform they must be 
identicaL Thus, if we find a continuous function y(t) 
from a given y{p\ for example by the processes described 
above which use the Table of transforms, this is the only 
continuous function which has y{p) for transform. 

It follows also from Lerch’s Theorem that, if is a 
constant, and 

y{p) = 0, p> k, then y(t) = o, o . (2.21) 

provided y{t) is continuous ; and if y(t) is not continuous 
it can only differ from zero at isolated points. 

•3. SoTne additional theorems 

Theorem III. If y{p) is the Laplace transform * of 
yif), and y{t) is continuous | and tends to y^ as t^ o, then 
i the Laplace transform of dyjdt is py(p) — yo, (3,1) 

If in addition, dyjdt is continuous and tends to as t~^o> 
then 

d^y 

the Laplace transform of ^ is p^y{p\ — py^ — 3;^. (3.;^) 

And if d^yjdt^ . . . dyjdt^ are continuous and dyjdtr y^, 
^ 3> * • ^ as t — >- 0, then 


the Laplace transform of 




“ P^yiP) - P^-^y, - - -pyn-^- yn-x- ( 3 - 3 ) 

These results follow immediately by integration bv* 
parts. To prove (3.1) we have 

- 4 >- 

functions concerned have Laplace 
°f such as for which the 
^ f not converge however large p is chosen, and 
transforms. Such functions “e 
£ IsS T? rt® to which the theory 

by the of § 24, TWeL IX.' ohtfxined 

“ ^“lue, at certaiix 

points, extra terms arise correspondmg to these points. 
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since as o, and ye~^^’^o as 00, 

provided p is sufficiently large. 

Similarly, to prove (3.2), vre have 




Ay , 
Jo 


= -yi-pya+P^- 

Repeating the process n times gives (3.3). 

Theorem IV. If y{p) is the Laplace transform of 
y{t)y then 

I 

the transform of 1 y{f)dt^ is . (3.4) 

Jo p 


For 


rco ci r T *1°^ T foo 

I e^^Au y{f)df — e~^^\ y(t*)dt' — I e'~'‘^^y(t)dt 

Jo Jo* L P Jo Jo ^Jo 



since the integrated part vanishes at both limits. 

The results of Theorems III and IV can often be used 
to obtain new transforms : as a simple example it follows 
from (1.7) and Theorem IV that 

a . I 

^ , — 5r is the transform of 1 sin at dt = -(1 — cos at), 

p(J>^^a^) Jo o' ’ 

Theorem V. If y{p) is the Laplace transform of y{t)y 
and Cxi is a positive constant^ then 

— y ^ is the transform of y {cat), . (3.5) 

For I = — J e~^^l^^^'y{t')dt^ ^^y 

This simple result is useful when it is desired to arrange 
the algebra of a problem in dimensionless form. For an 
example of its use see § 14. 

4. The solution of ordinary linear differential equations 
with constant coefficients. 
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Writing for d^yjdf^^ the problem is to solve the 
differential equation 

. . .+Cn-J^y-{-c^y=f{t), t>o, (4.1) 
where are constants, f(t) is a given function 

of and Dy, . . are to take the given values 

yo^ yv * • •» yn-^i when t = o. 

We multiply (4.1) by ^"^Vand integrate with respect to 
t from o to 00. This gives 

^'^e-P*{Dy + CiD"-b' + .. .+cj}dt= fip), (4.3) 

where f(p) is the Laplace transform of f(t), which can be 
written down* from Table I and Theorem I for the 
common types of function which usually occur. 

On the left-hand side of (4.2) we use Theorem III, 
which gives ^ 

pny _ + • ■ . +:yn-i) 

+ Cip'^-^y — + . . . + J'n-a) 

+ • • • + c„-ipy — c„-.iyo + =/(/>)• 

Or (pn + + . . . + c„)y 

=f(P) + c„.iya + c„_^(pyo + 3 'i) + • • • 

+ +P^~^yi + . • . + yn-x)- ( 4 - 3 ) 

The equation (4.3) is called the subsidiary equation 
corresponding to the difierential equation (4.1) and the 
given initial conditions. It is formed from the differential 
equation (4.1) by replacing D''y, r = o, on the 

left-hand side, by py- on the right-hand side, f{t) is 
replaced by its transform /(p), and terms containing the 
initial conditions are added according to the rule : 
for Dy -we add 
for „ py^+yi 

for DSy „ +^>3'i +P2 

for D’-p „ p'-ip, + pr-y^ + . . .+ 

* The solution for a function whose Laplace transform is 
not known, or even which has no Laplace transform, can be 
obtained by the use of Theorem IX, § 34. 
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The subsidiary equation (4.3) gives y(p)y the transform 
of the solution y(i) of the problem. To find y(f) from 
f(p), the methods of § 2 are used. Since, by the unique- 
ness theorem of §2, only one continuous function y(t) 
corresponds to a given f(p), it follows that this will be the 
unique solution of the problem. 

Example i. (D^ + 1)3; = 0, o 
with y = yoy Djy = y^^, when ^ — 0. 

The subsidiary equation is 

(P^ + i)y=pyo +yv 


Thus 

And 

Example z. 
with 


_ pyo + yi 
y p^+i' 

= 3^0 COS ^ + yi sin t, 
(D“ + 4D + 8)j = I, t> o 
y — 0, Dy = I, when t = o. 


Since the transform of i is i/p, the subsidiary equation is 
(j ,2 4. 4^ 4. 8)5; = 1 + I. 

The function y corresponding to this has been found in 
§ 2, Ex. I, and for a similar but more general case in 
§ 2, Ex. s- 


Example 3. (D + 1)^3; = sin zt, t > o, 
with y and Dy zero when ^ = o. 

The subsidiary equation is 

The function y corresponding to this y has been found in 
§ a, Ex. 2. 

Example 4. (D^ + n^)y ^ cos nt^ ^ > o, 
with y ^ yo>^y ^ yi^ when i ^ o. 


2 
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The subsidiary equation is 


Thus 


y=: 


Using (i.n), (1.7), (1.8), 


, pyo + yi 

it follows that 


y = — sin + Vo cos 4 - sin nt. 
zn n 


Example 5. (D + i)y = t> o, 
with y = y^^ when i = 0, 

Using (1.5) and Theorem I, the transform of 
is found to be zl{p + i)^. Thus the subsidiary equation is 

(^ + i)j=~-^+ 3 ;o. 

i.e. y = — — ? L - -^0 

^ (p + iy^p + i 

Then, using (1.5) and Theorem I, we find 

y = + y^e-K 

5* Simultaneous ordinary linear differential equations with 
constant coefficients. 

Here, proceeding as in § 4, using the transformation 
procedure on each of the equations in turn, gives a 
system of algebraic equations, the subsidiary equations, 
for the tr^sforms of the solutions. Suppose, for example, 
it is required to solve a system of n second order equations 

2 + *r.D + =fr{t), 

r = 1, 2, . . n, t> o, (5.1) 
w^re the f^{t) are given functions of t, and and T>y„ 
^ •. • *> to take the values u^ and when t — o, 

Multiplymg the equations (5-i) by integrating 

With respect to t from o to 00, and using (3.1) and (3.3) 
gives the system of n subsidiary equations 
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2 + Cr,)y, =fr{p) 

»-i n 

+ 2 Ki(^“s + ®») + ^rjWs]. »' = I, . . ., «. (5.2) 

These have to be solved for y^, . . 57^, and then 

Vi, • • •> 3'n ^re determined as in § 2. 

Example. (D“ + 2)a: — Dv = i-) . 

Djc + (D2 H-2)3 , = o/'^> ° 

to be solved with x = Xq, = y = Dj = o, when ^ = o. 
The subsidiary equations are 

{p^ + 2)® —py = J -\-px^ 

px + [p^ + 2):y = «(,. 

Solving we find 

r. ■ + l) + 2 I , pi^Xg — 1) ^p{ 2 ,X(f—l) 

p{p'^ + + 4) ^P 3 (P^ + 1)'^ 6 (p^+ 4) 

^ ~ (i>^ + iW + 4) 3 \p^ + I ~P~+^) ' 

Thus » = ■! + 1(2*0 — i) cos t + ^(2*0 — i) cos 2f 
y = j(2*o — i) sin t — -^(2*0 — i) sin 2t. 

6. Comparison with the ordinary method of solving differential 
equations^ and with Heaviside's operational calculus. 

The ordinary method of solving the linear differential 
equation of order n, 

,i>{D)y = (D" + + . . . + c„)y =f(t), (6.1) 

consists of finding the complementary function 

2 . . . (6.2) 

where are arbitrary constants, and the 

. . .j are the roots of (/>(p) — o (provided these are all 
different, if they are not the complementary function 
takes a slightly different form provided for by the theory) 
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and adding to this a particular integral determined from 
^(D) and f(t) by certain rules. This gives the general 
solution of (6.i). To find the solution satisfying given 
initial conditions, it is necessary to solve n algebraic 
equations for the constants A^, . . ,, A^. 

In the Laplace transformation method, the subsidiary 
equation (4.3) is 


y = “I" involving the initial conditions}. (6. 3) 


To determine jy from (6.3) we have to put the right-band 
side into partial fractions,* and for this the roots 

of (j){p) = o are needed ; thus these must be calculated 
in either method, but when they have been found, the 
work by the Laplace transformation method is much 
shorter ; for example, if are all different, 

y can be written down by (2.8), while with the older method 
of solution the n algebraic equations for the constants 
Aj, , . of (6.2) still have to be solved. Thus the 

advantage of the Laplace transform method over the 
ordinary one increases enormously as the order of the 
equation to be solved increases. 

The same remark applies to simultaneous linear differ- 
ential equations — the greater the number and order of 
these, the greater is the advantage of the Laplace trans- 
formation method, _ 

The part of y{t) derived from the term f{p)l 4 ^{p) of 
(6.3) is a particular integral of (6.1), but it may differ from 
the particular integral found by the rules of the older 
theory by terms of the complementary function. 

To see the connection between the present method 
and Heaviside’s operational method, we again discuss 
(6.1). Heaviside usually considered only the case in which 
the initial conditions were zero, and f{t) was his “ unit 
fimction ”, which is zero for t <io and unity for f >• o, 


* If is not one of the simple functions dealt with hitherto, 
it may not be possible to put /(p)/<^(p) into partial fractions, but 
the solution can always be found as in § 24. 
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and is represented by the symbol 1 . On the left-hand side 
of (6.1), D is replaced by^ ; this is regarded as an operator 
(not as a number as in the Laplace transformation method). 
Thus (6.1) leads to 

<l>{p)y = 1 


or 



(6.4) 


which is the operational expression for the solution y. 
The solution itself is then derived from (64) by Heaviside’s 
expansion theorem which is similar in type to (2.8). 

For this problem, the Laplace transform of the solution is 



and this differs from (64) only by the extra factor p in 
the denominator. In all cases the Heaviside operational 
expression is just p times the corresponding Laplace 
transform. 

As remarked in §1, the Laplace transformation pro- 
cedure is designed for, and applies immediately to, the 
‘‘ initial value problems ” of dynamics and electric circuit 
theory, in which differential equations have to be solved 
for times t > 0 with given initial conditions at ^ = 0. 
In other fields, such as the theory of elasticity, “ boundary 
value problems ” arise, in which a differential equation 
has to be solved in a region, say o <x < I, and the 
solution has to satisfy some conditions at the boundary 
a! = 0, and some at the boundary x 1 . The older 
method of solution, as the sum of a complementary 
function (6.3) and a particular integral, applies equally 
well to both types, since, in a boundary value problem, 
there will still be n conditions, some at :« == o, and some 
at X — I, to determine the A^, . . ., A„. The Laplace 
transformation can still be used, since the initial values 
y„ may be regarded as arbitrary constants ; some 
of these may be fixed by the conditions at == o, and the 
remainder can be determined from the conditions at x = Z 
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As an example, we consider the deflection * of a uniform, 
uniformly loaded, beam of length Z, freely hinged at the 
origin, and clamped horizontally at the other end, x = Z. 

. Writing D for differentiation with respect to x, the differ- 
ential equation for the deflection is 

EID^y = 0 <x <. 1 ^ . . (6.6) 

where E, I, and w are constants. This has to be solved 
with X = = 0 at a; = o, and x = Da? = o at a? = Z* 

We solve (6.6) with values o, o, oi y^ V>y, T>^yt 
respectively, at A!=.o, regarding y ^ and 3^3 as unknown, 
constants to be determined from the conditions at a? == Z. 
Writing y for the Laplace transform of y with respect 
to the present independent variable a?, the subsidiary 
equation is ^ 

p^y^^+pbi+y^^ 


Thus 


ZVX^ , VoX^ , 


The conditions y = Dy = o at * = Z, give 
wP y^P 

zoP 


Thus 




wP 


48EI’ 


y 3 


_ 3^1 

8EI’ 


and y = - sZ*^ + Px). 

7. The justification of the solutions of §§ 4, 5 from the 
pure-mathematical point of view. 

In deriving the subsidiary equations of §§ 4 and 5 some 
assumptions were made, which, from the point of view 


* This method is well suited to problems of this type, par- 
ticularly with concentrated or variable loads. For other examples 
see § 21 ; Jaeger, Math. Gaisiette, 23 (1939), 62 ; Pipes, Journ 
App. Phys.j 14 (1943), 486. 
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of Strict pure mathematics, would have to be defended. 
For example, it was assumed that y and its derivatives 
had Laplace transforms. To make the solution quite 
rigorous further discussion would be needed, but this 
can be avoided by checking the results obtained by the 
methods above, that is to say verifying that they do satisfy 
the differential equations and initial conditions of their 
problems. This can be done quite generally, and since 
the process itself is purely algebraic^ and adds nothing 
to the theory of the subject it is sufficient to state the 
results here. 

For the problem of § 4, the ordinary linear differential 
equation of order n with constant coefficients, it can be 
verified * that the solution obtained by the Laplace 
transformation method does satisfy the differential equation 
and initial conditions. 

For the problem of § 5, simultaneous ordinary linear 
differential equations, it is possible to give a similar general 
verification that the solutions obtained by the Laplace 
transformation method satisfy the differential equations, 
and that they satisfy the initial conditions, provided that 
the determinant of the coefficients of the highest powers of 
D in the terms does not vanish. This restriction is re- 
quired, not because of a flaw in the method of solution, 
but because of a fault in the specification of the problem. 
To see what it implies, consider the system of equations 
(5.1) : in this case the determinant in question is 

^nl • • • ^nn 

If this vanishes it is possible, by adding together appropriate 
constant multiples of the equations (5.1), to form at least 
one equation containing no second derivatives. Thus, 
instead of having a system of n second order equations, 
which could be solved with 2 n initial conditions, we have 
in fact a system of n — i (or less) second order equations, 

* C. and J., §§ 34, 35 ; Doetsch, loc. cit., Kap. 18. 
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and some first order or algebraic equations, and for such 
a system only zn — i, or less, arbitrary initial conditions 
can be prescribed. If the problem is stated with 
initial conditions, there must be relations between these 
or there will not be a solution. 

Problems of this type occasionally arise in practical 
applications, also similar difiiculties occur in switching 
problems in circuit theory. These questions are discussed 
in § 20 . 

EXAMPLES ON CHAPTER I 


I. Verify the following Laplace transforms : 


yiP) 


1 

<3 

1 


p 

(p — ^)(p “ 


I 

P(P^ + a^) 

' 

I 

' p^ip^ + a^) 

— sin at) 

I 

a6(&“ asin&O 

(P^ + a^)(p^ + b^) 

P 

(P^ + a'^){p^ + b^) 

rr— — sfcos at — cos bt) 
b^ — a^ 

I 

I 

Pl(p + a)» + 

«‘+ J3' j5V(a'‘ + /3*) 


sin \fit + arg(a + 


Solve the following differential equations with initial values 
^ 0 * yu . . of y, Dy, etc. 


(i) (D -h fi)y I, ^ + (^0 “ 

(ii) (D® + 5D + 6)y = i. 

b = -1 + (33^0 -h yi - + (i - 3^1 - 2yo)e-^°^3- 
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(iii) (D® + n^)y = sin a)£. 

r o) + 3;j(do^ — M^) . , 

^ • 2^ ■ sm nt + 3^0 cos 


i sin a)£ 


25 

]■ 


n{m^ — n^) 

(iv) (D + i)^3^ =: sin t, 

[y = (i + yo + yi)ie-* + (-i + ya)e~‘ — i cos t]. 


3. If 

with 


d^yjdx^ ~ Xj o < X < I, 
y == dyjdx = o at a; = o and x = I, 

y = + 2/® — 3x1^). 


show that 
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series with a leaky condenser of capacity C and leakage 
conductance i /G, the subsidiary equations are 

G + Q> 

. . ( 8 . 10 ) 

Returning to (8*6), we remark that it involves the 
initial current in the inductance and the initial charge 
on the condenser (these are related to the initial kinetic 
and potential energies, respectively, of the current). 
Further, since the terms involving them appear in the 
right-hand side added to the transform V of the applied 
voltage, it follows that the effect of initial charge Q on 
the condenser is the same as that of an applied voltage 
whose transform is (*— Q/Cp), that is of a constant voltage 
(— Q/C) ; similarly, the effect of initial current 1 in the 
inductance is the same as that of an applied voltage whose 
transform is Lt, this will be found in § 19 to be an irn- 
pulsive voltage Ll8(i). The same remark is true in 
general — ^that initial conditions may be treated as voltage 
sources of the appropriate types. 

9. Examples on L, R, C circuits. 

Example i. Constant voltage, E, applied at t ^ 0 to 
L, R, C in series, with zero initial current and charge. 
Here the subsidiary equation (8.6) is 

+ R + M 

Thus, 

T — ® ® r 

+ {RIL)p + (i/LC))- mp+t^y+nr ^ 

where _ R , _ i , _ i 

^ “ 2L’ ” “ LC ^ “ LC 4L2 

The notation (9.3) will frequently be used later. 

* For the exceptional case L — R = o and related problems, 
see § zo. 


+ R + -- ? + Ll 


Q ■ C(I+ 0 ) 
G + C^ 
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From (9.2), using Theorem I and (1.7), 
(1,9), respectively, we find 


I 


jg 

= sin nt. 

nL 


if > 0 


( 1 - 5 ). 


E 4.17 

rZjLj 


if = 0 

if == — < o. 
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Example 2. A condenser of capacity C charged to 
voltage E is discharged at t ^ o through an inductive re- 
sistance Ly R, 

Here Q = CE, I = o, V = o, and (8.6) becomes 

(Li) + R + ^)l= -|, . . (94) 

which is the same as (9.1) except for the negative sign. 
This implies that the current starts flowing away from 
the positive side of the condenser. 

Using (8.5) and (94), the transform of the charge on 
the condenser is found to be 

-__CE CE ' _ CE(p + 2ix) 

p p(LCp^ + RCp + x) + 

where the notation (9.3) has been used. Then, in the 
case > o, we have 

Q — CE^'^^l^cos ^ sin 


Example 3. Steady current E[Ris flowing from a battery 
of voltage y Ey through a parallel combination of condenser 
C and inductive resistance Ly R (Fig. 3), when at t — o the 
switch S is opened. It is required to find the potential across 
the condenser terminals. 

After the switch S has been opened, the circuit consists 
of L, R, C in series. The initial charge ^ on the condenser 
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is CE, and the initial current 1 in the inductance is — E/R, 
the minus sign being chosen, since, from the conditions 
of the problem, this is flowing away itom the 
side of the condenser. With these initial values (o.o) 
becomes 


(Lf + R + ^)I 


R p 


Thus . 


E(j> + 

WP + 


where the notation (9.3) has been used. Hence, con- 
sidering only the case > o, 

I = — . 5 e-/^^(^cos nt + ^ sin nt^ . 


s 



Fig. 3. 


If = Q/C is the voltage across the condenser ter- 
minals, it follows from (8.5) that 

. Q E E(Li> + R) 

® ~ C ~ ^ Rp(LCp^ + RCp + i) 

Er^ + (R/L) - (i/RC)] 
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This is the simplest type of circuit breaker restriking 
problem in which it is required to find the voltage across 
the contacts of a switch after they have been opened. 
In this case the voltage across the switch S is E — 
Another method of solving this problem is given in Ex. 2 
at the end of this Chapter. 


Example 4. Alternating voltage E sin (cot applied 
at if = o to an L, R, C circuit with ^ero initial current 
and charge. 

Since sin {(nt + 0) = sin wt cos (j> + cos cot sin cf), its 
transform is 


p sin <56 + ^0 cos <56 
+ 0)2 


(9-S) 


Thus, using the notation (9.3), the subsidiary equation 
(8.6) gives in this case 


Y _ E sin ^ cos 

^ L * + <^^)[{p + 


We consider only the most interesting case n®> 0. 
Here to find I we may either put (9.6) into partial fractions 
with real quadratic denominators, or notice that its 
denominator is 


(p — ia,){p + im){p + — in){p + M + *«). 

and use (2.9). The latter method gives * 

T Er— cu^sini^i + «a)®cos(ji ^ 

LI. Ztu>[(lJ, + twp + n^] J a j 

E f (- + inY sm <j> + wj- fj. + in) cos 

■^Ll 2m[(- u + 2«)2 4- cu2] 


-f Conjugate) (9.7) 


* Notice that if only the steady state part of I is required, only 
the terms of (a. 9) arising from the zeros ± iw of the denominator 
of (9.6) need be written down. The other zeros — ti± in give 
transient terms. See also § 13. 
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To reduce this to real form, ‘ we use (^.14) as iti §2, 
Ex. 5. Thus, writing 

^[2(jLco + i{oj^ — — n^)'\ = . (9’^) 

6t) 

(jLL^ + — n^) + zifin = 

[(/x^ ~ n^) sin 0 — (ji>o) cos <^] + cos cf) — 2/xw sin 

(9,7) becomes 

F EZ 

I == I sm (a)^ + ^ - 8) + sin + B,), ( 9 - 9 ) 

where Z and B may be written out explicitly as the 

and argument of the complex number on the left of (9-^)i 

and similarly Z^, Zg, Sg. 

10. Electrical networks 

It is convenient to regard a complicated electrical net- 
work as built up of a number of circuit elements of th< 


^5 Is 




type studied in § 8, with their terminals connected in som( 
specified manner. Let L^, R^, Q be the inductance 
resistance, and capacity in the r-th of these elements 
(Fig. 4), let V,. be the voltage drop over it, let I,, be th€ 
current in the element, with initial value I,., and let Q 
be the charge on the condenser, with initial value <!i,. 
To be quite general mutual inductance must be incluLcied 
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SO let be the mutual inductance * between the r-th 
and j-th elements. 

Then the equations for the r-th circuit element are 

+ + ^ + v„ . (lo.l) 

(. 0 . 2 ) 

Proceeding as in § 8, the subsidiary equation for the 
r-th circuit element is found to be 

+ LL + 2 M„l,- (10.3) 

fi+r e^r ^rP 

where a:, = + ^, . . (10.4) 

.... (10.5) 

From the form of the network, further equations can 
be written down connecting the and in different 
circuit elements. Firstly we have Kirchhoff's First 
Law ] 

I = o, at a junction. . , (10.6) 

This leads to 

I = o, at a junction. , . (10.7) 

Kirchhoff’s Second Law states that, for any closed circuit 
in the network, the sum of the externally applied voltages, 
V, in the closed circuit, equals the sum of the voltage 
drops, Vy, over the circuit elements included in the closed 
circuit. That is, writing ^ to denote summation round 

0 

Mf, will have a sign depending on the way in which the 
coils are wound, and on the convention of sign we set up for 
the Currents in the different circuit elements. 

t Usually labour can be saved by Specifying the currents, 
eitheri directly or by the use of mesh currents, so that they satisfy 
(10.6) automatically. 

3 
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a closed circuit 2 V, 

0 0 

and thus • ‘ ‘ (io.8) 

0 0 

Using (io.8) in conjunction with (10.3) we have 

o'- a=ff 0 0 ^ a+r ^ 

2 , (10,9) 

By choosing suitable closed circuits, sufficient equations 
of types (10.9) and (10.7) can be found to give all the 1^, 
"\^en voltage V, applied' to the terminals of any net- 
work with zero initial conditions, produces at the terminals 
current I given by 

:^(/>)I = ?, . . . (lo.io) 

we call z the “ generalized impedance of the network 

viewed from these terminals ; * for example, defined 
in (10.4) is the generalized impedance of an Rr> 
circuit. Now for a single circuit element of a network 
with no mutual inductance and with zero initial conditions, 
(10.3) becomes 

zX=%- ■ • • (lO.Il) 


That is, I,., and are connected in the same way as 
current, voltage, and resistance in direct current theory. 
Thus in this case the generalized iinpedance of the network 
can be written down by combining the of the circuit 
elements which comprise it according to the ordinary laws 
for combining resistances in direct current theory. For 
e;xample, the generalized imped^ce of a combination of 
n generalized impedances 5^2, . * in parallel is 

given by 


I 

z 


= T + r + - • 


(10.12) 


* z(ia)) is the ordinary complex impedance of alternating 
current theory. 
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Again, the generalized impedance of the ladder network, 
Fig. 5, can, by repeated application of the result for two 
resistances in parallel, be written down as the continued 
fraction 



Fig. s- 


It follows in the same way from the remarks at the 
end of § 8, that the current in any branch of a network, 
due to initial current 1 in an inductance and initial charge Q 
on a condenser in branch, is given by 

z{p)l^hi-{tllCp) . . ( 10 . 14 ) 

where z{p) is the generalized impedance of the network 
viewed from a pair of terminals in the branch. 

n. Examples on electrical networks 

Example i. The two circuits of Fig. 6 are coupled by 
mutual inductance M, Voltage V{t) is applied to the 
primary at t = o with zero initial conditions. It is re- 
quired to find the secondary current {the primary resistance 
is neglected). 

Here the subsidiary equations (10.3) give 

+ M/.I, = 

+ (h2p + 1^2 + ~ 
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Mp7 




' (L^La - M^)p^ + R^L,p + {LJC^)' 


(ii.i) 


V 


r 


3 


o ^ 

M I L; 

11 


IfV^ 


and 


Fig. 6. 

: E, constant, so that ^ = 'Ejp, this becomes 


Is = 


I« = 


ME 


{L,L,-M^)p^ + R,L,p+(LJC,y 

ME 


(ii.a) 


where a = 


^(LiL,-MT 

^2^1 02 
V p ■ 


-<xt 


sin 


La 


2(LiL2-M2)’ C2(LiL2- M2) 

The case = M^ is discussed in § 20. 

Example 2. In the circuit of Fig. 7 the condensers are both 
charged to voltage E, when, att = o, the switch S is closed. 


L 


B 






L G 

— I 


i i-i, 
c 


D 

Fig. 7. 


H 


We choose the currents as shown in Fig. 7 to satisfy 
Kirchhoff’s first law. Then the subsidiary equation 
(10.9) for the closed circuits ABD and BGHD are, 
respectively, ^ 
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Solving for I, we find 
__ E ^ ^ 


_ E f I + Vs 

where w 2 =i/LC. 

■ Therefore 


1-V5 

y+i(3+V5K/’ 


I- -=V(^){y5^»”«[l(3-Vs)l‘ 

Example 3. /tz the circuit of Fig, 8 it is required to find 
the current in the inductance L caused by prescribed 
current I supplied at the terminals AB, starting af f = o 
•with zero initial conditions. 



From (10.12) it follows that the generalized impedance 
of the combination of L, R, C in parallel is 
I _ 

I I ■“ RLj!) 

+ R + RLCp 2 + + R. 

Also the generalized impedance of the inductance L i? 
Lp, Since the same potential is applied to both we hav 

RLC^)^ + Lp + R^' , * ■ 
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For example if I = i, for i > o, so that I = i/p, 
and, if > o, 

where a = <f> = tan-^ (^/a). 

Example 4 . Bridge networks. 

We consider a general bridge network represented by 
the scheme of Fig. 9, and suppose known voltage V to be 


■ i)(RLCp2 + L^> + R)’ 
sin (^^ + <;6), 



H 


applied at i = o with zero initial conditions. It is required 
to find I5, which is the galvanometer current in most 
applications. Choosing the currents as shown in Fig. 9 
to satisfy Kirchhoff’s first law, the subsidiary equations 
for the closed circuits ABDHA, ABFA, BDFB, re- 
spectively, are 

-t- — Is) = ^ 

+ ^ 616-^312 , =0 

Is) -^4(12 I5) ^sls ^ 


^lll 
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Solving for I5 ,we find * 

^ + -^2 o 

^ A 

^2 ^4 

= |(« 2 ^ 3 -%^ 4 ). 


I 

o 

o 


where 


A = 


+ ^2 


d 


^2 


% 

^2 


““^3 ^6 

*^4 % '^4 -^B 


(11.4) 


(II-S) 


The way in which the bridge can be used for measure- 
ments is determined by the form of in (11.4) 

as a function of p. If, by suitable choice of the circuit 
constants, this can be made to vanish for all p^ we have 
I5 == o, and thus, by (3.21), 

15 = 0, £>0, 


whatever the applied voltage V. As an example of this 
case consider the bridge in which consists of inductance 
L and resistance in series, 3*2 of resistance Rg, z^ of 
resistance R3, and z^ of capacity C and resistance R4 in 
parallel. Here 


Z-i — \jp -|- R^, ^2 — ^ 2 > -^3 — ^3? -^4 


Thus 


^2^3 — ^1^4 = 


R2R3(^ ~b R 4 QP) — 

I + R 4 CP 


R4 


I + ^/S^P 
R4(^.P + Ri) 


Then if R2R3 = RiR4> 

and R2R3C = L, 


we have I5 = 0 for f > 0, whatever the form of the applied 
voltage may be. 


^ It is hardly necessary to point out the advantage of using 
determinants, particularly when, as here, only one of the currents 
is needed. 
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As an example of a bridge whiqh cannot be balanced in 
this way, consider the case in which 

5^1 — -|“ Ri> *^2 ^ ^2) ^3 ~ 

and consists of resistance in series with a parallel 
combination of resistance rg and capacity C, so that 

I ^ _ h'^2p'P + + ^2 

+ r,Q. + i ■ 

and 5?2 ;s^3 — 

LCriy2j>^+jft{(Riri— RaR3)Cya+L(ri+ra)}+{Ri(^i+ra)--R-8l^a} 

Clearly (ii.6) cannot be made to vanish for all values of p 
by choice of the circuit constants. However, we may 
use the circuit as an alternating current bridge by requiring 
that the steady state part of I5 is to vanish when the applied 
voltage V is sinusoidal and of frequency wjzTT, Antic- 
ipating the results of § 13, this requires that the real and 
imaginary parts of the numerator of (11.6), with^ replaced 
by io)y shall both vanish, that is 

L(^i + '‘ 2 ) = C;ra(R2R3 - Vi) \ (rx.-f) 

and LCr^r^oj^ = ^^{r^ + r^ — R2R3/ * ^ 

Alternatively the bridge may be used as a ballistic 
bridge ; in this case a battery of voltage E is applied at 
t — o with zero initial conditions, so that V == E/p in 
(ir.4), and the requirement is that the total charge passed 
through the galvanometer is to be zero, that is 

{ 00 

\dt = o. 

0 

By Theorem VIII, § 22, this requires 

lim Ig := 0. . . . (11.8) 

>0 ^ 

It is easy to verify from (11.5) that A does not vanish 
when p = Oy so from (11.4) with V = E/p, and (11-6), 
it follows that for (i 1.8) to be satisfied we must have 


R2R3 ““ Ri(^i + ^2)1 
R,r2^C = L(r, + r2) J 
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In practice, instead of applying a battery to an initially 
dead circuit, it is usual to open the battery circuit when 
steady current is flowing from the battery. In that case 
the impedances of Fig. 9 must be precisely 

specified, and the circuit of Fig. 9 (with the arm AHD 
removed so that 1 ^ — — I^) has to be studied with known 
initial values of the currents in the inductances and the 
charges on the condensers. 


12. The natural frequencies and damping factors of a circuit 
The transforms of the currents, say Ii, . . I„, in the 

various branches of an electrical network are found by the 
solution of equations of types (10.9) and (10.7), and 
always take the form 


I r- 


n, 


(12.1) 


where 0 (p) is the same for all I^, it is in fact the determiuant 
of the coefficients of the in the equations and depends 
on the circuit only. The numerators fr{p) depend on r, 
and involve the initial conditions and the applied voltages. 

Before we can find the I^ we need the roots of the period 
equation ” 

0 (p) = o, . . . (1S.2) 


and we know that to a root — /x i zV of (12.2) there cor- 
responds a term sin (v£ + S) in the I,., that is, a term 
of period 27 r/v and damping factor 

The essential point is that the period equation (12.2), 
whose roots give the natural frequencies and damping 
factors of the circuit, is always found in the course of the 
calculation -of any of the I,., but the contain in addition 
the complete solution of the transient problem. Thus in 
§ ir, Ex. 2, it follows from the first expression for I that 
the natural frequencies of the circuit are w[|(3 ± V s)? I 
In demonstration problems, such as most of the exarnples 
solved here, *P(p) is usually arranged to have obvious linear 
or quadratic factors. In practice *^(p) may be a poly- 
nomial injp of degree n with coefiicients given numerically. 
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an example of such a problem is given in § 14. Here 
we give some references on the methods available for the 
solution of algebraic equations. 

For cubic and quartic equations algebraic methods are 
available, described in text books on algebra.* * * § For the 
cubic, some tables and charts are given in Jahuke-Emde.l 
Numerical methods can be used for cubic and quartic 
equations, and must be used for equations of higher degree. 
Real roots can be found from first principles, using a graph 
and improving the rough values obtained in this way, 
either by interpolation (an example of this is given in § 14) 
or by Newton’s method (W. & R., § 44). A cubic has 
always one real root which can be found in this way ; 
dividing by the corresponding factor leaves a quadratic. 
The more sophisticated methods for finding real roots are 
Homer’s- method (C. S., § 475, W..and R., § 53) and the 
root squaring method (W, & R., § 54), the latter is prefer- 
able. Problems on resistanceless circuits lead to equations 
in with real (negative) roots which are easily solved by 
the above methods. For equations having several pairs 
of complex roots, the modified root squaring method seems 
the most satisfactory ; J another method, and also a method 
of successive approximation to a complex root, is given by 
Frazer and Duncan. § 

If any of the roots of the period equation have positive 
real parts, the solution will contain terms which increase 
exponentially with the time ; in this case the system is 
said to be unstable. If all the roots have negative real 
parts the system is stable. Clearly it will be useful to 
have a result which will indicate whether a system is stable 

* Charles Smith, Treatise on Algebra^ Edn, 5 (1921), referred 
to as C. S. ; Turnbull, Theory of Equations (Oliver and Boyd, 
1939). Whittaker and Robinson, Calculus of Observations 
(Blackie, 1924), referred to as W. & R. 

t Tables of Functions (Edn. 2, Teubner, 1933). 

t Brodetsky and Smeal, Proc. Cambridge PhiL Soc., 22 (1924), 83. 
Note that when there are several pairs of complex roots the method 
given in W. & R. is unsatisfactory. 

§ Proc. Roy. Soc., A 125 (1929), 68. 
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or unstable, without the labour of solving the period equa- 
tion. One such criterion is Hurwite's * which states that 
the roots of the equation 

= o 

all have negative real parts if all the determinants 


% ^0 

> 



0 

, • * . 

% 

^0 

00 .. 

, 0 

^3 ^2 




% 


^3 


Oi ao 

0 



^6 




^2n- 

-1 ^2n- 

-2 

• 


1 

(12.3) 


are positive, where, in writmg down the determinants, 
is put equal to zero if r > w. 


13. Steady state alternating currents 
It was remarked m § iz that if voltage V is applied to 
any circuit at t = o, the transform of the current I,, in any 
branch of the circuit will always take th.6 form 


T -V /’■(^) I Sr(p) 


■ (13-1) 


where fr{p) and 0 {p) depend on the circuit only, and 
g^{p) involves the initial conditions. Now suppose 


so that 


V == E sin + ^,) 

^ CO cos 6 p sm (j> 

V — h j 5 y 

p^ + 


(13-2) 


and (13. i) becomes 

T p cos sin 4 )fr{p) , grip) 


(13-3) 


If the system is a dissipative one, the roots of = 0 
will all have negative real parts, and thus will give rise to 
transient terms. If only the steady state altern ating current 


* Frank-v. Mises, Differentialgleichungen der Physik (Edn. a, 
1930), Vol. I, p. 1 6a. This is equivalent to Routh's Rule, 
Advanced Rigid Dynamics (Macmillan, Edn. 4, 1884), Chap. vi. 
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is required, we need only evaluate the two terms of (^-9) 
corresponding to the zeros ih of These give 

^/- (co cos ,j> + tVsin ^ Conjugate} 

^ Z^) 0 — Sr) 

where — Z,(co)e*‘’r. 


Thus, if the transform of the current in any part of R 
circuit due to applied voltage V and any initial conditions 
has been found, the steady state alternating current can 
be deduced immediately. An example has been given in 
(11.7). As another we find the steady state alternating 
current in the secondary of Fig. 6 due to the primary 
voltage E sin (aJ^ + (f>). Here (ii.i) and (13.Z) give 
I _ MpE (o) cos ^ sin </>) 

' “■ (P^ + co2)[(L,L2 - M^)p^ + + (Li/C^] - 

The steady state alternating part of I2 is 
/ M2<iuE(a> cos <j> + ioj sin _• _ ^ 

\2ioy[-w\LJu^-M^)+(LJC^)+RjL.^icS\ + Conjugate 





McoE 


cos (^cot -j— 


5 ), 


where Z = +[(LJC^) _ (L^L^ - 

and 8 = arg {[{LJC,) - - M^)] + iR,L,c^y . 

If z(p) is the generalized impedance of a dissipative 
circuit looking in at a pair of terminals, and voltage 
is applied at these terminals with, zero initial conditio las, 
the transform of the inpilt current will be 

E 

{p - iio)z(py ■ 

and the steady state current will be 
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14. Numerical treatment of more complicated circuits 
In the problems previously considered, the denominators 
of the transforms of the solutions have had obvious linear 
or quadratic factors, so that explicit expressions for the 
solutions could be written down immediately. In more 
complicated circuits, the denominators may contain ex- 
pressions of the third or higher degree in p which have 
no simple factors. If the circuit constants are given 
numerically, the zeros of the denominator can be found 
by one of the methods mentioned in § 12, and the solution 
completed ; if the circuit constants are unspecified, families 
of curves can be obtained which will show the behaviour 
of the circuit under any circumstances. 


±L 


I 


c 



Fig. 10, 


In this section we study, from both these pomts of 
view, the current I in the circuit of Fig. lo, due to unit 
voltage applied at f = o with zero initial conditions. In 
the usual way we find 

, + 1) . (X4.1) 

^ “ RCLiCiP“ + Li(^ + ^ 


As an example in wHch v?o-"S 

suppose these are . 


Lx = o-i, henrys ; 
becomes 

1 = 


Cx= io-», C = 2-S X 10 
R = 10,000, ohms. Then 

io-^(p* -h loio) 


, 15 * 


p3 1-4 X ^ 

To evaluate I, the toota of the ^ 

p 3 _j_ 1.4 X 10®^® + lo’-V + O 4 X 


farads ; 

(14.1) 

(14.2) 

(14-3) 
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are needed. Putting p = icPx in this, gives the more 
manageable equation 

f{x) = + i- 4 »^ + A? + 0*4 = o. . (144) 

This cubic can be solved by any of the methods of § 12, 
but it is very simple to work entirely from first principles. 
From a graph it is found that there is a real root at 
approximately = — 0*77. This root can be found 
more precisely as follows : we find /(— 0*77) = 0-00353, 
/( — 0-78) == — 0-00279, between 

— 0*77 and — 0*78 ; further, if f(x) ran linearly between 
its v^ues at these two points it would vanish at — 0-776. 
By repeating this process, as many places of decimals as 
desired can be found. Working with three places of 
decimals only, we divide f(x) by (x + 0-776), and obtain 


f{x) = (« -f o‘776)(jc^ + 0-624A? + 0*516) 

= h + 0’776 )[{x + 0*312)^ + (0-647)^]. (i 4 ‘S) 

Returning to (14-2) we find, u^sing Theorem II as in 
§2, Ex. 5, that the function whose Laplace transform is 


IS 


<2 + a® 






+ 


V- 


( a ^^) 2+^2 

where 8 =arg(f 2 -f a~ 7 j 2 + 2 ffT 7 )+arg(a— 


sin {y]t — 8), ( 14.6) 

(14*7) 


Using the values a = 0-776 X lo^ ^ = 0*312 X 10®, 
r) = 0-647 X 10^, a ~ 10^®, found above, and writing 
i = io~®T, so that T is the time in micro-seconds, (14.2) 
and (14.6) give finally 

ioT==:2*527e-o-0776T 4, 1.^33^-0*03121 gin (0-0647T-- 1-485) 

(14.8) 

To plot I from (14.8) we notice that o*o647T = tt/S 
when T = 6*07, so that steps of 6*07 in T correspond to 
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increases of in the argument of the sine, and only 
eight values of the sine need to be looked up. The cal- 
culation may be set out as follows, A and B being the first 
and second terms on the right hand side of (14.8), 


T 

q*o776T 


A 

0*03I2T 

6 -it 

Tji— d' 

sin 

B 

10*1 

0 

0 

1 

2*527 

0 

1 

— 85*i« 

—0*996 

-1*527 

1*000 

6 *o 7 

0*471 

0*624 

1*577 

0*189 

0*828 

~62*6“ 

-0*888 

— 1*127 

0*450 

12*14 

0*942 

0*390 

0*986 

0*378 

0*685 

-*40*1° 

—0*644 

—0*676 

0*3x0 

i8*21 

1*413 

0*243 

0*614 

0*567 

0*567 

— 17*6° 

—0*302 

—0*263 

o* 35 X 

24*28 

1*884 

0*152 

0*384 

0*756 

0*470 

4 * 9 ® 

0*085 ‘ 

o*o6i 

0*445 

30*35 

2*355 

0*095 

0*240 

0*945 

0*389 

27 * 4 ® 

0*460 

0*274 

0*5x4 

36*42 

2*826 

0*059 

0*149 

1*134 

0*322 

49 ’ 8 '’ 

0*765 

0*378 

0 - 5^7 


It very frequently happens in practice that it is not the 
transient behaviour of a circuit for given numerical values 
of the circuit constants which is needed, but a general 
description of the transient behaviour for any values of 
the circuit constants. The circuit of Fig. 10 will now 
be studied from this point of view. The first thing to 
do is to express (14.1) in a form involving only simple 
and fundamental quantities, such as dimensionless ratios, 
time constants of parts of the circuit, etc. Thus, writing 


fiA.i) becomes 

r,T_ + ^ . (14.XO) 

+ nbp^ + 

A further simplification comes from expressing (i 4 -io) ^ 
terms of A = P /«i. so that the solution is obtamed m terras 
of tijt by Theorem V. Thus writing 

p^=plni, k = nlth, . ■ (H-ii) 

, TjT ^ . . (14.12) 

we have - p^s kbpj_^ + Pi + k 

The right hand side of (14.12) oontams the two. ^ime^ 
sionless, parameters 6 and ; to study the general behaviour 
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of the circuit we can regard one of these as fixed at some 
^bitrary value, and see how the solution changes as the 
other takes all possible values. We shall give h fixed values, 
and for each of these allow k to run from o to oo, that is, 
say, R to decrease from oo to zero. 

First we have to study the position of the roots of the 
period equation 

+ khx"^ -|- re + ^ = o . . (14*13) 

for constant h and k >0. We could give k a number of 
values and solve (14.13) for each of these, but this would 
be very laborious, and also we could not be sure that the 
values of k chosen would fairly represent all types of be- 
haviour in the circuit. Both these difficulties can be 
avoided by the following procedure : the cubic (14.13) 
has always one negative real root, call this — a, and 
suppose that it has also a pair of complex roots — ^ ± ^17 
of modulus r. Then we know that * the sum of the roots 
of (14.13) equals — kh, the sum of their products two at 
a time equals i, and the product of the roots equals — K 
Or, written in ffill, 

. . (I 4 -I 4 ) 

+ 2 a^ = I . . . (14-^5) 

ar 2 = A , . . . (14^16) 

where + 77^. . (i4-i7) 

If we regard a as a parameter, and solve (14.14) to (14.16) 
for and k in terms of it, we find 




i+a^ a(6 — i) _a(i+o^^) 


(14.18) 


That is, — a, — ^ — ^^), with the values of r® 

and ^ given by (14.18), are the roots of (14.13) with the 


* Charles Smith, Treatise on Algebra. If a^, . . are the 
roots of 

+ . . . -f a^-xx 4- an ~ o 

— ' a^f 2!oC;lCC 2 — fl2> ^ ““ ^^3, ■ • 

ttitta . . . a„ == 


then 
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value of k given by (14* i8). If the value of ^ given by 
(14.18) turns out to be greater tlian. r^, the same argument 
shows that there are three real roots, 

- a, and - I ± _ ^2) 

for the corresponding value of k. 



Fig. II. 


In Fig. II the position of the complex root ~ f + J**? 
of (14.13) is plotted as a function of the parameter a, 
for the values 2, 4, 8, 9, 10 of h. When a = 0 this root is 
i for all values of h, and as oc increases the curves shown 
are traced out. The positions of the roots for the values 
0*2, 0*4, 0-6, i-o, 2-0 of a are marked by successive crosses 
on the curves, k is given by (14.18), and, when there is 

4 
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a complex root, increases as a increases, that is, as we 
move round the curves. If ^ 9 there is always a pair 

of complex roots of (14.13). lih — ^ there is a real triple 
root for h ^ 1/(3 Vs) at oc = i/Vs J for all other values 
of k there is a pair of complex roots. If & > 9 there are 
three real roots for some values of k ; these are indicated 
by plotting (dotted in Fig. ii) half the sum of the real 
roots other than a as abscissa, and half their difference 
as ordinate, measured below the axis. 

^ diagram such as Fig. ii gives an immediate 
indication of the way in which the roots of (14*13) vary 



with k and h. To study the way in which the current I 
derived from (14.1a) varies with k and we select typical 
values of h from Fig. ii, for example, the transition value 
6 = 9, a value of b less than 9, and a value greater than 9. 
On each of these curves we choose a few representative 
points md calculate I for these using (14.6), 

In Fig. 12 the resulting set of curves for the case b — q 
IS shown. Here RI is plotted against n^t for the values 
0-05, o-i, 0-2, i/V3» 4 j 8, 16, 32 of a, to which 
correspond the values 0-049, ^*152, 0-192, 0*^70, 

0-469 0-901, 1*784, 3*558 of k. Clearly these curves, 
together with two or three such sets for other values oi 
^ fairly complete indication of the way in 
which the transient behaviour of the circuit depends on 
k and h. 
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15. Circuits which include vacuum tubes 

Transient problems in valve circuits are easily treated 
by the Laplace transformation method, provided, of course, 
that the valves are operated in linear parts of their char- 
acteristics. Only triodes are considered here, but similar 
considerations apply to other valves. 

The characteristic equation of the triode supplies a 
linear relation between the anode current, 1^, and the 
anode, grid, and cathode voltages V^, V^, V,,, respectively, 
namely pla = (V„ - V,) + - V„), . (iS-i) 

where p is the anode slope resistance, and jit the amplifica- 
tion factor. In transient problems, the quantities V^, V^, 



Vo, la appearing in (15.1) are the variable parts of the 
voltages and current caused by the transient behaviour 
being studied ; they are superposed on the constant vol- 
tages and anode current caused by the anode, and grid bias, 
batteries. The grid current is negligible in normal cir- 
cumstances, and here is always assumed to be zero. In 
any problem it must be decided whether the interelectrode 
capacities in the valve itself can be neglected— if not, 
they are easily included as shown below. By virtue of 
(15. i) it is possible to replace the valve by an equivalent 
circuit element consisting of resistance p and a voltage 
inserted between anode and cathode; 
this is the common practice, but it is a little more funda- 
mental to work directly from (15.1) as will be done here. 

As a first example we consider the type of problem that 
occurs in the amplification of transients : in the circuit 
of Fig. 13 interelectrode capacities are neglected, and 
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steady conditions prevail due to grid bias and plate batteries 
(not shown in the figure), when at = o a known addi- 
tional voltage is applied to the grid ; it is^ required to 
find the change, I„, in the anode current. Since Yg and 
the changes in currents and voltages it causes are super- 
posed on the steady values, we do not need to know the 
latter, and need only write down equations for the changes, 
all of which vanish at t = o. All the currents and voltages 
marked in Figs. 13-16 are to be understood in this sense. 
We notice that (15. i) implies a convention of sign, namely, 
if > Vo, is flowing in the direction of the arrow in 
Fig. 13 ; the same current also flows from the cathode. 



The subsidiary equations for the impedances ^0 

. (15.2) 

S’ola = • • • (IS'3) 

the signs being determined by the convention stated above. 
Also the transform of (15,1) is 

plfl = -- (i + ' (^S-4) 

Solving (15.2), (15.3) and (15.4) we find 

{p + + (l + l^)^o}‘^a — * (^S-S) 

and from this can be found if is known. 

As a more definite problem we take the case of resistance- 
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capacity coupling, Fig. 14, and find the value of Vg (which 
would be the grid voltage of a second valve) corresponding 
to constant voltage E applied at i — o to the grid of the 
first valve. Here we have, in the notation of Fig. 13, 


^0 

Also 


G + Cp’ “ 1 + C2(Rx H- Ra)/ 


"^2 — F.2I2 — 


RiRgCajpIa 
I + Ca(Ri + R2)^ 



jtxE(G-|-’ Qp)RiR2C2 

[I p(G+ C^] [i + CaCRi-f ^2)p] + Ri( I + R2Ci^)(G + CpY 

(iS-6) 



where (15.5) has been used in the last line of (15.6). Since 
the denominator of (15.6) is a quadratic in p, it is easy to 
write down V2. 

To illustrate the effect of inter-electrode capacity and 
a rather more complicated system, we consider the circuit 
of Fig. 15. Here grid-plate capacity is included in ^3, 
grid-cathode capacity in and anode-cathode capacity 
in ;^2- We seek the anode voltage due to voltage V 
applied at i = o at the terminals AB when steady con- 
ditions prevail. 
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Choosing currents as shown in Fig. 15 to give zero grid 


current, the subsidiary equations are 

= . ( is -7) 

^(^1 + ^2 ” = ^{7 * * (^5‘8) 

^2X2= . . (15.10) 


and the transform of (15.1), which is, in this case, 



Fig. 16. 


Solving we find 

y 'V^^2(p (I'i'ia) 

whence can be found if the impedances z are known. 

Finally we treat an oscillator circuit (Fig. 16) from the 
present point of view. Here there is a closed system in 
which the grid voltage is not prescribed from outside, but 
is free to adjust itself. As always, there are steady currents 
and voltages due to plate and grid bias batteries, and the 
1, 1^, are departures from these : clearly one solution 

is to have all of these zero, and the problem is to find how 
this solution behaves if the system is disturbed. To fix 
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ideas * suppose there is a change of current t in the 
inductance L 2X t — o, then the subsidiary equations are 

(L^ + R)I = ■?, + Ll . (15.13) 

. (15.14) 

(Li^ + Ri)(L + I)= . (15.15) 

• (lS-l6) 

The equations (15.13) to (15.16) can be solved for 
h Ia» ^a> ^ terms of t and the circuit constants. All 

the solutions will have for denominator 

A(^) =:Lj:> + R o o — I 

— I o Cp — Cp 

hj^p hip + Ri I o 

O p — I •— /M, 

=LL,C(i+i^)^^+C^^[p(L+Li)+(i+/x)(RL,+LRi^^ 

+jP[p^(^i + R) + Li + CRiR( I +p)] +p + Ri, (15*^7) 

and the solution contains terms in where a^, r=i, 2, 3, 

are the zeros of (15.17). If all the have negative real 
parts the disturbance dies out gradually and the system is 
stable ; if any of the has a positive real part the system 
is theoretically Unstable, actually it executes an oscillation 
of amplitude determined by the limits to the linearity of 
the characteristic. The condition for this, which is the 
condition for the circuit to act as an oscillation generator, 
is obtained by applying Hurwitz^s criterion, (12.3), to 
(15.17); it is 

C[p(L + Li) + (i + /^)(RLi + LRi)] 

[pC(R+Ri)+Li +CRRi(i +i^)] < CLL,( i +p){p+R^). 
16. Filter circuits 

The essential feature of such circuits is that a number of 
precisely similar circuit elements are arranged to form a 

This is done solely as an illustration, all that is needed is 
the determinant A which does not involve the initial disturbance. 
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repeating pattern. Mathematically this implies that we 
can write down a general relation, called a difference 
equation, connecting the transforms of the currents in 
successive elements, and by solving this we get a general 
expression for the current at any point. 

As a definite example which includes many cases of 
practical interest we consider the circuit of Fig. 17. Here 
there are m impedances and w — i impedances ^ , 
arranged as shown. Voltage V is applied at the beginning 
^ough impedance and the whole is terminated by 
impedance jstq. We require the current at any point due 
to the voltage V applied at ^ = o with zero initial conditions.* 



Fig. 17. 




The subsidiary equations for the successive meshes of 
Fig. 17 are 

(2sr+s:')l,+i+a:r+2=o. r = 0, i, . . (16.2) 

+ (^0 + • • • (^^'3) 


The equatiops (16.2) which connect the transforms of 
any three consecutive currents are a system of difference 
equations. We seek a solution of them of type 
where A is independent of r. Substituting this in (16.2) 
gives 

z — ( 2 Z + == 0, 

or + e-'^) — (zz + z') — 0, 

z' 

that is cosh 0=1+--. . . (16.4) 

zz 

* If the initial conditions are not zero, as in the discharge of 
a^ charged filter circuit, there mil be terms on the right hand 
sides of the subsidiary equations which depend on the particular 
impedances z and z" in the circuit. 
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. Thus, if e is defined by (16.4), a solution of (16.2), 
in fact the general solution * of it, is 

Ir — Ae’’® + 

It follows from (16,4) that sinh^ ^9 === and, since 

both signs of 6 occur in (16.5)? the square root with the 
positive sign may be taken, which gives 


sinh ^0 = ► (i6-6) 



Fig. 18 (a). Fig. 18 (6). 


To find the constants A and B in (16.5), this must be 
substituted in (16.1) and (16.3), which give two equations 
for A and B, namely 

~ ze^)A + ^ — ze''^yB = 

A(;2^o + ^ + B(;sro + == o* 

Solving these and substituting in (16.5) gives finally 

(g+gp) sinh (m — r)6—x smh (Tn—r— i)^ ^ 

(g-f-gj)(g-f gp) sinh sinh sinh (m 

(16.7) 


This general expression simplifies greatly wh.en_ 
and Zi have values corresponding to the common tenpina- 
tions. The network of Fig. 17 may be regarded as built up 
of either T elements (Fig. iSa) or X] elements (Fig. 180) 


* The theory of difference equations is analogous to that of 
differential equations. (16.5) corresponds to the comploment^y 
function, and, if there had been functions of r on the right 
side of (i6.»), a particular integral would have had to be i^luciecl. 
A brief reference to the theory is given m Piaggio, Dtnerential 

Equations, p. The Laplace transformation method be 

applied to difference equations, cf. Gardner and Barnes, ioc. ctt. 
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connected in tandem, and the general expression (16.7) 
takes specially simple forms when circuits built up in this 
way are open or short circuited. For example : 

© ^0= — o : Fig, 17 with no series impedance, 

short circuited at the far end, 

I == V cosh {m-r^ i)e 
" 2^ sinh (m - i)d sinh ^ 0 ' ' 

(ii) ^0= CO, Zi — 0. The same, but open circuit at 

e far end, ^ ^ 

^ 2^ cosh {m - sinh ^< 9 ’ * 

(iii) A circuit of m, T sections, with 
no additional series impedance and short circuited at the 
far end. Here, using (16,4) we get 

f cosh (m — r )6 , ^ n 

^ 2: sinh ^ sinh ‘ * 

(iv) z^ = Zq— CO. The same as (iii), but with 
open circuit at the far end, 

T _ "V sinh {m — r)d 
^ z smh 6 cosh md' * * 

(v) Zi = z\ Zq ^ z' 2,z. A circuit of m + i, ff 
sections, with no additional series impedance, and open 
circuited at the far end, 

I = ^ sinh {m~r + i)d 

2Z cosh (m + i)d sinh ^ 6 ' * i • / 

If ar, S'', Zq, Zi are known functions of p, (16.4) gives 
cosh 0 as a function of p, and using the formulae 

cosh nd = 2^-1 f cosh d — cos ( cosh 6 — cos 

\ 27 z/ \ znj 

. . . (^cosh e - cos ~ (16.13) 

sinh nd = 2«-i sinh 9 (^cosh 9 — cos ^cosh 9 — cos — ) 

. . . ^cosh 9 — cos — (16.14) 


n 
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it appears that (16.7) can always be expressed as a quotient 
of polynomials in p. Both (16.13) and (16.14) 
any positive integer w, even or odd. It is often advisable 
to use them to see the exact form of the denominators of 
the transforms of the currents, though actual evaluations 
of the currents 1^ are best carried out by working in terms 
of 6 as in the examples below. 

Example i. Constant voltage E applied to a circuit 
of tn, T sections in which z' = R, z = i/Cpy with no 
additional series impedance, and open circuit at the far 
end. 

Here (16.4) gives 

cosh 0=1 + ^RC^, . . (16.15) 


And, using (16. ii), 

» _ CE sinh {m — r)d 
’’ ~ sinh 6 cosh md 


(16.16) 


To evaluate we need to know the zeros of the de- 
nominator of (16.16). These can be found from (16.13) 
with the value (16.15) of cosh 0, or by noticing that cosh mO 
is zero when 


w — I. (16.17) 


“ Taking more values of than those in (16.17) would 
result in counting the same zeros over again. By (16.14), 
sinh 0 is a common factor of the numerator and de- 
nominator of (16.16) so 0 = 0 need not be considered 
in using Theorem IL 

The values of p corresponding to the values of 0 given 
by (16.17) using (16.15), 



^(Z5+I)7r-| 


., TW—I. (16.18) 


To evaluate I,, it is best to use the form (2.8) of Theorem 
II, and for this we need 
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rd 

\jdp ^ 




= [^(smhecosh«0).|]^_^^^ 

_ fiRC ^ + 7w sinh ^ sinh 

L^ sinh ^ oc. 


= JRCtw sinh 


sinh ^ 

(2^ + i)z7r 


= |^zw2RC(— )«. . . . , , (16.19) 

Then, appl3dng (2.8) to (16.16), using (16.19), gives 


m-:. 




(m - r){2s + i),r 


sm 


|«“0 


2 ,m 


(16.:^:^) 


WzR ^ 27 « ^ ^ 

Example 2. Voltage V = sin oj^ is applied at jf “ 0 
to a filter of T sections, in which z' = Lj&, z = ilCjpy 
with open circuit at the far end. 

Here (16.4) gives 

cosh ^ = I + ^‘LCp^, . . {16,-2, i) 

And from (16. 1 1) 

j __ Cpo) sinh {m — r)Q 
^ '{p^ + (^'^) sinh d cosh mQ* 

The zeros of the denominator of (16.22) are + ico, 
and the values oi p corresponding to 

a (25 -|- l)z 7 r 

2^^ ^ “ I- (^^-^ 3 ) 

Using (16.21) these are seen to be 

P = ±i^s> . .... (16.24) 

A = . (r6..s) 

and it is assumed that none of these equals co. 
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Proceeding as in Example i we find 

r~ (sinh d cosh mQ)\ == LCm/S^. 

\_ap 

And thus from (2.8) and (16.22) 

( 2 ^ + i)r7r 
2 m 


(16.26) 


2ct) 

ml. ^ 






j8S 

Co) sinh {m 


cos cos ■ 
— r)0(fa)) 


Usinhe(fa,)cosh« » W " + Conjugate} 
where 0(zcu) is the value of B corresponding to ^ = 
The second term of (16.27) steady state 

period 27 r/co, To evaluate it, we see from (16.21) 
cosh 6 {ico) == I — 

and this gives rise to two possibilities : 

(i) if 0)2 < 4/LC, e{ia>) == zS, 

where 


(16.27) 


ICO, 


part of 

that 

(16.28) 


8 = cos”^ (i — -ILCcu^), 


and in this case the part of of period 27r/aj is 


Coj 


sin {m — r) S 


cos ojt. 


sin 8 cos mB 
(ii) if 0)2 > 4/LC, d(ico) == zV + A, 

where A = cosh“^ (^LCco^ ~ i), 

and in this case the part of 1^ of period Z'lrjco is 


sinh A cosh wA 


cos cot. 


(16.29) 


(16.30) 


As r is increased, (16.30) decreases exponentially, while 
on the other hand (16.29) oscillates. That is, if co^ > 4/LC 
the term of period 27 r/co dies out exponentially as we move 
along the circuit, while ii co^ C. 4/LC this term oscillates 
trigonometrically. The circuit is, of course, a low pass 
filter, and we have found both its transient and steady state 
behaviour for applied voltage sin cot. 


17. Vibrations of mechanical systems 
Problems on the motion of bodies of concentrated inertia, 
connected by shafts or springs of negligible inertia, lead 
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very readily to systems of ordinary differential equations 
which have to be solved with given initial conditions. In 
order to make these differential equations linear, the dis- 
sipative forces are restricted to the case of resistance to 
motion proportional to the velocity, though in some 
problems Coulomb friction (constant and opposite to the 
direction of mption) can be treated. Rotating systems 
only will be considered here, since they are of perhaps 
a little more interest than the corresponding problems in 
linear motion, and lead to precisely similar equations. 

We write throughout this section, J for the moment of 
inertia of a wheel about its axis, 9 for its angular displace- 
ment from a standard position, D for differentiation with 
respect to time, ^ ^ ^ 

for angular velocity, kco for resistance to motion proportional 
to angular velocity, T for torque, and A for the stiffness of 
a shaft, that is the torque required to produce unit relative 
angular displacement between its ends. 

The equation of motion of a wheel of moment of inertia 
J, acted on by torque T, and with resistance to motion 

(JD + A)a. = T. . . (17-2) 

Also the torque T required to produce a relative dis- 
placement 6 between the ends of a shaft of stiffness A is 

T = AR . . . (17.3) 

Using these equations we can write down differential 
equations for the motion of any combination of wheels 
connected by shafts of known stiffness, A number of 
examples of the differential equations which arise in such 
problems is given below : clearly these can be solved just 
as in § 5 with given initial values of the co and 9 , 

Example i. A shaft of stiffness A has one end fixed, 
and on the other a wheel of moment of inertia J (Fig. 19(a)), 
Torque T, a given function of time, is applied to the wheel, 
and there is resistance kco to its motion. 

If 9 is the displacement of the wheel from its equilibrium 
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position, the retarding torque due to twist of the shaft is A^, 
and the differential equation for the motion of the wheel is 


CID+k)c=^T-Xd. . . (17.4) 



Fig. 19 (a). Fig. 19 (b). 


Equations (17.1) and (17.4) can be solved with given 
initial values of 0 and in the usual way. If preferred, 
o) can be eliminated ; this gives the second order equation 


(jOa + + A)0 - T. . . (17.5) 



Fig. 20 (b). 


Example z. Three wheels are connected by shafts of 
* stiffness A and Ai as shown in Fig. 20(0), in which the 
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symbols used all have the same meanings as before, By B^, B2 
being displacements of the wheels from a reference position 
in which the shafts are unstrained. Torque T is applied 
to the wheel J. 

Here the differential equations for the motions of the 
wheels J, J^, respectively, are 

(JI)+k)co^T + X(B^^e) . . (17.6) 

(J,D + k,)co^ = - Xie, -^ 6 )+ X^( 6 , - d,) (17.7) 

(J2^ ^2)^2 “ \(^2 ' ^l)> • • (^7*8) 

and CO = D^, <0^ = = D^g. , . (i7*9) 

This system can be solved in the usual way with given 
initial values of 6 , B^y d^y a>, co^^ Clekrly the problem 
of any number of wheels connected by shafts of given 
stiffnesses can be solved in the same way ; if all the wheels 
have the same moment of inertia and the same resistance 
to motion, and all the shafts the same stiffness, the sub- 
sidiary equations form a chain of the type met in § 16 
and are treated in the same way. 

Example 3. A geared system. 

First it is necessary to study the type of complication 
introduced into the system by 
gearing. Suppose a gear of 
radius drives a gear of radius 
7*2, their angular velocities being 
and 0)2 in the directions of 
the arrows in Fig. 21. Then 
if and B^ are the angular dis- 
placements of the two gears 
from their original positions 

»'A = »'A • • (17-10) 

rjcoj = rjcoa. . . (17.11) 

At the point of contact of 
the wheels there is a driving 
force P on the wheel r^y and 
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an equal, oppositely directed, reaction force P on the 
wheel r^. Then if T2 is the driving torque on the wheel rg, 
and Ti is the reaction torque on the wheel 

T 2 = ^2?, Ti = r^P 

and thus r{Y-^ = r^Tg. . . . (17.12) 



Fig. {a). 



(17.11) and (17.12) form a pair of linear connections 
between the angular velocities and torques in the two 
shafts. More complicated systems, such as differential 
gears, give rise to linear connections in the same way. 

As an example we consider the system of Fig. 22 {a) 
in which torque T is applied to the wheel J which drives 
a gear of radius and moment of inertia through a 

5 
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shaft of stiffness A. The gear meshes with a gear of 
radius and moment of inertia J^, which drives a wheel 
of moment of inertia Jg through a shaft of stiffness 
With and Tg defined as in (17.12) the equations of 
motion are 

(JD+*)to = T + A(0i-0) . (17.13) 

(J,D + - «) - . (17.14) 

(JaD + == Ai(^3 — ^2) + Tg . (17*^5) 

(J3D + — ^2)- • (17-16) 

These equations, together with (17.1), (17.10), (17.11), 
(17.12) are sufficient to give the motion for any initial 
conditions consistent * with (17.10) and (17,11). 

In many cases it is possible to set up an electric circuit 
problem which corresponds f to a given mechanical 
vibration problem, in the sense that both lead to precisely 
the same systems of differential equations. This does not 
make the problem any easier to solve, but it gives an inter- 
esting physical insight into it, also the corresponding 
electrical system may be set up and studied in the laboratory, 
and finally, since the theory of electric circuits is more 
highly developed than that of mechanical systems, it may 
be useful in suggesting further developments in that 
field. 

To see the sort of correspondence J which exists between 
mechanical and circuit problems, compare equations 
(17.2), (17.3), (17.1), respectively, from which the differ- 
ential equations for mechanical problems were built up, 
with the fundamental equations from which the differ- 

* Initial conditions inconsistent with these correspond to 
forcing into mesh at £ o a pair of gears inconsistently moving 
or displaced. 

t For a large number of such equivalent problems, cf. Olson, 
Dynamical analogies (van Nostrand, 1943). 

X The correspondence described here is not the only one ; 
thus, instead of (17.17) we might have used the equation for 
a leaky condenser, namely (CD + G)V — I, and in this way 
another, equally useful, system arises ; cf. Miles, Jfourn. AcousU 
Soc. America, 14 (1943), 183. 
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ential equations of circuit problems are built up, namely : 

(LD + R)I = V, . . (17.17) 

for an inductive resistance with applied voltage, V ; 

g.Q = V, . . (17.18) 

for a condenser with applied voltage, V ; and 

DQ = I, . . (17.19) 

for the connection between charge and current. It appears 
that (17.2), (17.3), (17.1), and (17.17) to (17.19) are the 
same sets of equations in terms of different quantities, 
the correspondence between these being set out in the 
Table below, in which the corresponding quantities for 
linear motion of a mass are added for completeness. 

Electrical. Rotational. Linear Motion. 

Charge, Q. Displacement, d. Displacement. 

Current, I. Angular velocity, co. Linear velocity. 

Inductance, L. Moment of inertia, J. Mass. 

Resistance, R. Damping coefficient, k. Damping coefficient. 
Capacity, C. Reciprocal stiffness, 1/ A. Reciprocal stiffness. 
Voltage, V. Torque, T. Force. 

Using this correspondence we can now set up electric 
circuits which have the same differential equations as 
the mechanical examples considered earlier. For example 
the circuit of Fig. 19 {h) with inductance J, resistance A, 
and capacity i/A, with current co, charge B on the con- 
denser, and applied voltage T, gives the same differential 
equation (17.4) as the mechanical problem of Fig, 19 [a). 

Similarly Fig. 20 {b) is the equivalent circuit * to Fig. 
20 {a), so that the differential equations found by applying 
KirchhofFs second law to the closed circuits AA^BiB, 
A^AgBaBi, and AgAgBgBg, are (17.6) to (17.8) respectively, 
where By B-^y 6 ^ are the amounts of charge which have 
passed through the inductances J, Jj, Jg at time t, so that 

* Leaky condensers in Fig. 20(b) correspond to resistance to 
motion proportional to the relative velocities of the wheels in 
Fig. 20 (a). 
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0 -^ 0 ^ and ~ ^2 the charges on the condensers 
I /A and i/Ai respectively. 

Finally we consider the equivalent of the geared system 
Fig. zz (a). The linear connections (17. ii) and (17. la), 


i.e. 



(17.20) 


become in the electrical case 



(17.21) 


which are the connections between currents and voltages 
in an ideal ^ transformer of secondary and primary 
turns. On this understanding the equivalent circuit to 
Fig. zz {a) is Fig. zz {h), and it is easy to check that these 
lead to the same differential equations. 


18. Servomechanisms 

The ideas involved | may be illustrated by considering 
a rotation control system in wliich a massive body of 
moment of inertia J has to be driven so that its displace- 
ment Oq follows as closely as possible the prescribed dis- 
placement Q^{t) of a pointer. To do this, voltage 9 ^ — 0 ^ 
is fed into an amplifier-motor system which applies torque 
k{ 9 i — ^0) to the body, where A is a constant. This torque 
tends always to bring the output displacement 6 ^ 
into agreement with the “ input displacement ” 9 ^, and is 
proportional to the “ error ^ To provide 

damping for the system, it is usual to “ feed-back ” into 
the amplifier a voltage depending on the output velocity 
D^o, or on the error velocity D^. The simplest arrange- 
ment is to mount a generator on the output shaft which 
provides voltage which is fed back to the amplifier 

♦ Cf. Starr, Electric Circuits and Wave Filters (Pitman, 1934), 
p, 128. 

t Cf. Hazen, “ Theory of servo-mechanisms ”, Journ. Franklin 
Inst.f 218 (1934), 279 ; Minorsky, " Control Problems ”, ibid., 
232 (1941), 451, 519. 
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with negative sign. The torque supplied by the motor is 
then ki^i — Oq — and the equation of motion of 

the body is = k{e^ - Oo - . (i8.i) 

If the motion is started at ^ = o from rest at zero dis- 
placement, the subsidiary equation is 

Qp^ + kk^ + k)e,=^ke,, . . (18.2) 

or ^0 = .2 _L iT 2! ' T; — 2’ * * (^^- 3 ) 

" k-^ojrp + 


where 


• (18.4) 


Suppose the input is given a uniform motion d^ = 
starting at ^ == o, so that 

6 i^Q,lp\ (i8.s) 

then (18.3) gives 

Oi I ^i^iP "b k^co^ i) TiS 6^ 

p^ p p'^ hyOJ^p 0 }^ 

It follows that 


cos nt + — i) sin ni|, (18.7) 

where — oj^ ^ 

It is seen from (18.7) that the error Q ^ Oq — d^in this 
case is composed of a constant lag, — kiQ , , together with 
transient terms. 

If instead of feedback the feedback is arranged 

to be proportional to the error velocity, the equation 

of motion is ]D% = k{— 6 — kj)d) . . (18.8) 

or JD 2<9 + kkj)e + k 9 = - 

If the motion is started from rest at zero displacement, 
the subsidiary equation is 

(p2 k^oj^p + w ^)6 = — p^Oiy , (18.9) 
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where oy^ = k/]. If, as before, 6 i ~ Q,tj this gives 

s 

and ^ sin . . (18.10) 

n ^ 

In this case the error (18.10) consists of transient terms 
only and there is no permanent lag as in (18.7). 

In practice the voltage or is usually not 

fed back directly to the, amplifier, but applied to some 
circuit, and voltage is tapped off from some part of this 
circuit and fed to the amplifier. In this way rather more 
complicated transforms arise. 

Another complication of practical importance is the 
presence of time lag in some part of the system. To 
illustrate the effect of this suppose that, in the first case 
considered above, the voltage — ^0 — instead of 

being fed directly to the amplifier, is applied to an induc- 
tance L and resistance R in series, and that the voltage 
across the resistance is fed to the amplifier as before. 
Here there is a time delay determined by the time constant 
of the L, R system. The subsidiary equation for the 
current I in the L, R circuit, with zero initial conditions is 

(L^ + R)I = - (i + h^p%, . (18.11) 

and the voltage V across the resistance R is given by 

V=RL . . . (18.1a) 

Also the subsidiary equation for the motion of the body 
of moment of inertia J, with applied torque is 

. . . (18.13) 

From (18.11) to (18.13) it follows that 

Q ^ /^o y .X 

0 {p^ ^ <^%p ^ i)P> ^ 

where to® = kj], co^ = R/L. . . (18.15) 

The error 0 ^ may be evaluated, for given 6 ^,■ when the 
roots of the period equation 

{p^ + o)'^){p + coi) + oj\kyOJ-^ — i)p — o (18.16) 
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are knowA. Hurwitz's criterion (12.3) for the roots of 
this equation to have negative real parts, that is for the 
system to be stable, requires 

• • • (18.17) 

Thus if there is time delay in the system it is necessary 
for the amount of feedback to be greater than the minimum 
amount (18.17) system will not be stable. The 

way in which the roots of (18,16) and the values of 
vary as is increased may conveniently be studied by the 
method of § 14. 

19. The impulsive function 

In dynamics the idea of a very large force acting for 
a very short time and such that the time integral of the 
force is finite is extremely useful, and the “impulsive 
motion ” of dynamical systems is extensively studied. 
In circuit theory the ideas of an impulsive voltage (that is, 
a very great voltage applied for a very short time and such 
that its time integral is finite) or an impulsive current, 
are equally fruitful.* 

It is convenient to have a symbolic representation of 
such a function and this may be obtained as follows. 

Consider the function f(t) defined *|* by 

f{t) = 0, t <0 1 

= l/£, 0 <t<€l . (19.I) 

= 0 i > e J 

* For an interesting account of impulsive functions and their 
applications see van der Pol, Journ. Inst, Elect. Eng., 81 (1937), 
381. 

t Many other functions have the same properties, e.g. the 
triangular function, Fig. 23 (&), or the function 
(i /«V 5 f) exp (— tvs’*), 

Fig. Z3 (c), except that the latter has its maximum at t — o. 
It is more usual to choose a function symmetrical about « = o, 
such as Fig. 23 (c), but there is some advantage in the present 
applications in being able to think of, say, an applied voltage 
which is zero at £ = o and becomes very large immediately after 
this instant. 
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where e is very smail. This is very large in a very small 
region to the right of the origin (Fig. 23 {a)) and zero 



elsewhere, and, no matter how small e is, 


( 8 
8 

II 

M 

(19.2) 

Also, if /(^) is a continuous function 


fco T ra-he 

1 fit)cl,(t — a)dt = -\ f{t)dt 

J —00 f J a 


= /(a) + 0(.), . 

(19-3) 

where 0(6) is written for a term involving 6 as a 
and thus which tends to zero as 6 0. 

We now define the “ impulsive ” function h{t) to 
the properties of as e — >• 0 that is : 

factor, 

possess 

S(^) = 0 for # < 0, and for i > 0, except for a*^ 
vanishingly small region to the right of ^ = 0 y 

(194) 

TOO 

1 = I, 

J — 00 

(I- 9 - 5 ) 

J — 00 

(19.6) 


if/(i) is a continuous function. 
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It follows from (19.6) that 

[ =1, . . . (19.7) 

J 0 

that is, the Laplace transform of the 3 function is unity. 

We may now regard a blow of impulse P as a force 
P8(^) ; an impulsive voltage whose time integral is E, 
as a voltage ES(^) ; and an impulsive current which transfers 
charge Q as a current QS(i), and proceed to solve problems 
involving them in the usual way : the validity of the 
solutions so obtained is discussed later. 

Example i. A mass M attached to a spring of stiffness 
is set in motion at i == o from rest in its equilibrium 
position by a blow of impulse P. The equation of 
motion is 

MB^x + Moi^x = PS(i), 
and, using (19.7), the subsidiary equation is 

Thus X = (P/Ma>) sin cjot, . . (i9‘8) 

Example 2. An impulsive voltage EB(t) is applied at 
^ = o to an L, R, C circuit with zero initial conditions. 
Here the subsidiary equation (8.6) is, using (19.7)^ and the 
notation (9.3), ^ 

^ = L[(p + + n^- 

Thus, in the case > o, 

I = cos nt — fx sin nt). . (i9*9) 

Lw 

The solutions of Examples i and 2 above have been 
obtained by the ordinary Laplace transformation procedure, 
using (19.7) for the transform of 8(^), and, since any use 
of the 8 function must be regarded as symbolical only,* 

* The S function was defined as the limit of some function as 
€ -> o. Any mathematical operations on this function, such as 
differentiation or integration, would involve the interchange of 
order of two limiting processes and be very difficult to justify. 
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we must inquire whether these solutions are mathematically 
justifiable, • The simplest way to do this is to verify that 
they do in fact satisfy the differential equations and initial 
conditions of their problems. Consider the solution 

(19.8) from this point of view : it satisfies the differential 

equation ^£>2^ ^ ^ . (15.10) 

for the motion of the mass M for r > o ; also as ^ 0 

in (19.8), a: o, and Dx ->• P/M. Now P/M is the velocity 
with which, according to the theory of impulsive motion 
in dynamics, the mass M would be started by a blow of 
impulse P. Thus the solution obtained by the formal 
Laplace transformation procedure applied to Example i 
does in fact give the complete solution of the dynamical 
problem correctly, and in addition it avoids the necessity 
of evaluating the initial velocity due to the blow. It can 
easily be verified * that the same is true for a general 
dynamical system of n degrees of freedom set in motion 
by blows. 

Precisely the same is true for the electrical case in 
Example 2. Here, as a the value of I given by 

(19.9) tends to E/L, and we need an electrical analogue 
of the theory of impulsive motion in dynamics to show 
that an impulsive voltage ES(^), applied to an L, R, C 
circuit in which no current is flowing, will instantaneously 
produce a current of this magnitude. This can be deduced 
from the differential equation for current and charge in 
the circuit, namely 

(LD + R)I + g = ES(f). . . (I 9 -II) 

Proceeding as in dynamical theory, we integrate this 
over a very small time, r : the integral of 8(^) over this 
time is i ; that of the finite quantities I and Q is very 
small ; and that of DI is the difference between the final 

* By an easy extension of the method in Carslaw and Jaeger, 
loc. cit.y § 35. Jaeger, PhiL Mag. (7), 36 (194s), gives a 
complete verification for all the types of problem arising in this 
section and § 20. 
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value, ly, and the initial value (zero) of the current I. 
Thus LI/ = E, . . . (19.12) 

and the theory of impulsive motion gives the current 
ly = E/L, to which the solution obtained by the Laplace 
transformation procedure was found to tend as i o. 
It is easy to verify that (19.9) satisfies the differential 
equation of the problem for ^ > o, so it does in fact give 
the correct solution of the physical problem. 

It can be proved in general that for a complicated circuit 
with various impulsive applied voltages, the solutions 
obtained by the Laplace transformation procedure, treating 
8(i) as having (19.7) for transform, do satisfy the equations 
of circuit theory, and, as ^ o, the currents tend to the 
values which would be found from the equations of 
impulsive motion. 

zo. Jmtification of the solutions of electrical network 
problems 

It was remarked in § 7 that, for a system of linear 
differential equations, a general verification could be given 
that the solutions obtained by the Laplace transformation 
method satisfied the differential equations and initial 
conditions, provided a certain determinant did not vanish. 
The problems of circuit theory often give rise to a system 
of differential equations to which this verification applies 
immediately. In general, however, they give a combina- 
tion of differential or algebraic equations (10. i) and alge- 
braic equations (10.6) ; it may happen that the specified 
initial values of the currents and charges do not satisfy 
some of these equations, and in that case there will be an 
impulsive readjustment of these initial values. It is in 
fact true * in all cases that the solutions obtained by the 
Laplace transformation method, used as in § 10, giving 
each inductance its own initial current, and each con- 
denser its own initial charge, in the subsidiary equations, 
do in fact give solutions of the equations of the problems 

Jaeger, loc, cit. 
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for ^ > o, which as ^ o tend to the values which exist 
in the circuit after the impulsive readjustment of initial 
conditions, if this takes place. The general theory will 
not be given here, but only examples illustrating the 
different cases which can arise. 

Example i. Steady current E/R is flowing in the 
circuit of Fig. 24 with the switch S closed. At t = o 
the switch is opened. 


s 



Since the initial current is E/R in L, and zero in 
and Lg, the subsidiary equations are 


Solving we find 
E 


(Lp + R)I + L2i)l2 = ^-h|^ 

I = Ix + I.- j 


. (20.1) 


I 


LiL„E 


R{(LLi + LL„ + LiL,)/) -I- R(Li -|- L„)} 

And thus 


(20*2) 




L1L2E 


^ ^-R(Li+L2)f/(JXi'l-inL2+LL..) (20 'l) 

R R(LLi+LL2+LiL2) 

A . \ ^ H~ LL.) ^ ^ 

> O, > ^ ^ (20.4) 
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and, since this does not equal E/R, it appears that there 
has been an impulsive readjustment of current between 
the three inductances just after the switch was opened. 

To investigate what happens in this case, and to verify 
this result, consider the equations of the system for 3? > o, 
namely 

(LD + R)I Hh L^DL, - E, . . (20.5) 

L0DI2 — LjlDIi — o, . . (20.6) 

I2 + Ii = I. . . (20.7) 

The values E/R, o, o, which the currents I, I^, Ig had 
before the switch S was opened, do not satisfy (20.7). 
Suppose that these change to E^), and in a very 
short time r following the opening of S at ^ = o, then, 
integrating (20.5) and (20.6) from ^ = o to t=r, and 
neglecting the integrals of the finite quantities E, I, I^, Ig 
over this very small time, we find 

L{I(o) - (E/R)} + = 0, . (20.8) 

= 0. . (20.9) 

Also I<^), and satisfy (20.7), that is 

1(0) = 1^(0) 1^(0). . , jo) 

Solving (20.8) to (20.10) gives the new values of the 
currents, and, in particular, I(o) is found to have the value 
(20.4). If the voltage drops across the inductances were 
evaluated in the usual way, they would be found to contain 
S functions, that is, the finite changes of current may be 
regarded as being produced by impulsive voltages. 

Example 2. The condenser C in the circuit of Fig. 25 
(p. 78) is charged to potential E, and the condenser Cj is un- 
charged, when at i = o the switch S is closed. 

Flere the subsidiary equations are 

Cp p 


. (20.11) 
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Solving we find 

. _ EC(RCi/> + I) 

R(C-fCi)i. + i 

_ ECCi EC* I 

" ' C + Cl R(C + C,)* > + [i/R(C + Cl)]- 

TLu. I _ - ^S(,) - 

R 



Fig, 25. 


The term [— ECCi/(C + QL)]S(i) in (20.12) represents 
an impulsive current which instantaneously exchanges 
charge ECCl/(C + Ci) between the two condensers ; 
the second term represents the subsequent discharge of 
the condensers through the resistance. Whenever a 
network includes a closed circuit containing condensers 
only, there will be an instantaneous exchange of charge in 
this way between the condensers if these are incompatibly 
charged. If we had used the subsidiary equations for 
charge instead of those for current, this finite change of 
charge would have been found (in the same way as the 
finite change of current in Example i) but as we have 
worked with current, the impulsive current which transfers 
this amount of charge has appeared. The result (20.12) 
can be verified by writing down the differential equations 
for the charges on the condensers and proceeding as in 
Example i. 
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Example 3. Steady current E/R^ is flowing in the 
primary of a transformer, Fig. 26, and the secondary 



current is zero, when at i = o the switch S is opened. 
Writing Ij and for the primary and secondary currents, 
the subsidiary equation for the secondary is 

+ Ra)!^ = ME/Ri . (20.13) 

Also Ii = o, for t>o, so Ij^ = o. 

Thus from (20.13) 

lii = . . (30.14) 

Here as i o, Ig ME/R^Lg. There was a sudden 
change in the currents when the switch was opened 
which could have been studied as in Example i. 


Example 4. The problem of § ii, Example i, for the 
case of LjLg — (perfect coupling). . Here (11,2) becomes 
- ME 

In 


R2L1P + (L,/C,)’ 


and thus 


I2 == — --—eH/RaOj 


(20.15) 

(20.16) 


It is known from general theory that ; the 

limiting case LiL2 = is of interest for two reasons : 
firstly it may be very v nearly approached in practice ; 
secondly, if = M^, the degree of the denominator 
in I2 is reduced by one, here a quadratic has become linear 
and similarly a cubic becomes a quadratic, and the algebra 
is greatly simplified. 
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The significance of perfect coupling from the theoretical 
point of view * is best seen by considering the differential 
equations of the circuit of Fig. z6 : these are 


(LiD + Ri)Ii + MDI2 - E| 
MDIi + (L2D + R2)l2 = o/ 


(20.17) 


If L1L2 = multiplying the first of (20.17) by L2, 
the second by M, and subtracting, we obtain 

LgRJi — MRaIg = LaE. . (20.18) 

Thus in this case we can form an algebraic equation by 
elimination from the given differential equations, and if 
the given initial currents do not satisfy (20.18) they will 
change suddenly to hew values which can be calculated 
as in Example i. 

The examples above illustrate the types of problem in 
which there may be sudden changes of charges or currents 
on closing switches, and the way in which these could 
have been found has been sketched briefly. A general 
discussion along these lines can be given, which shows 
that the Laplace transformation procedure, applied in the 
usual way, does in fact always give the correct solution of 
the physical problem. 


EXAMPLES ON CHAPTER II 


1 . Alternating voltage E sin (nt + e) is applied at t = o to 
inductance L and capacity C in series, = i/LC. If the 
initial current and charge are zero, show, using (i.ii) and (i,.i2), 
that the current at any time is 

(E/2n'L){nt sin (nt +€)-}- sin c sin nt}. 

2. Constant current (E/R) is supplied for i > o to a parallel 
combination of inductive resistance L, R and capacity C with 
zero initial conditions. Show that the voltage drop across the 
condenser is 

E - Ee-^* {cos nt + [(R/2nL) - (i/«RC)] sin nt], 

* This is an example of the exceptional case referred to in 
§ 7, since the determinant discussed there is LiLg — for 
the system of equations (20.17). 
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where /a and n are defined in (9,3). Deduce the solution of 
§9, Ex. 3. 

3., In the problem of § ii, Ex. 2, show that the charge on the 
first condenser is 


(CE/2V5){(i + V 5 ) cos 7zi[K3 - Vs)]* . , ^ 

- (i ““ Vs) cos ; 2 t[i (3 + Vs)]*}* 


4. Regarding an imperfect condenser as capacity C in series 
with resistance Ri, the whole being shunted by conductance G, 
show that the subsidiary equation for voltage V applied to such 
a condenser in series with an inductive resistance L, R is 


{r^-^ ^ Q> + G(i + RiCp)}^ V + Li G(^+RiCi.)’ 

5. In the bridge of Fig. 9, consists of resistance Ri and 

capacity Ci in parallel, 2^3 consists of resistance R3 and capacity 
C3 in parallel, and js'a and 2^4 consist of resistances Rj and R4, 
respectively. Show that there is a balance for*all types of applied 
voltage if Q/Ca = R./Ra = Ra/Ri- 

6. In the circuit of Fig. 6, L1L2 = M^, and steady-state con- 
ditions prevail due to alternating voltage of peak value E and 
frequency co/ztt in the primary. At the instant of zero secondary 
current the condenser is instantaneously discharged, show that 
at time t after this instant the secondary current is 


ECaVLa 


V{Lx(i + co^Ra^Ca^)} 


o) sm cut 


RaC: 


■e-tlRiO,y 


7. In the resistance- capacity coupled amplifier of Figs. 13, 
14, with Zo — o, show that if the grid voltage is changed by 
E at £ = o, the change in grid voltage of the next valve is 

KiKa "h p(Ri + 1^2} 

8. A mass M hangs by a spring of stiffness Mw^. At t = o, 
when the mass is at rest in its equilibrium position, the other end 
of the spring is given a downwards motion a sin a>^. Show that 
the displacement of the mass from its equilibrium position at 
time t is 

-r sfo) sin nt — n sin a)i]. 


9. Constant Jvoltage E is applied at t ~ o to a condenser of 
capacity Ci in series with a leaky condenser of capacity C and 
leakage conductance G. Show that, if the initial charges are 
zero, the current is 


ECCi 

C + Q 


m + 


EGCi^ 

(C + c,) 


.e-G«/(04Ci). 


6 
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FURTHER THEOREMS AND THEIR 
APPLICATIONS 

21. Heaviside* s shifting theorem 
Theorem VI. If y{p) is the transform of y{t\ and 
a > o, then e~^'^y{p) is the transform of 

y{t — d)'R{t — a), . . 

where * H(0 = o, t<o\ , 

= 1, t>oJ ' • ^ 



Fig. 27 (a). 


Since H(^ — a) is zero if t < a and unity if t'> cl, the 
function y{t — <a:)H(i — a) is zero for i < a, and for t> a 
its graph is exactly that of y{t) for ^ > o, shifted to the 
right through a distance a, as in Fig. 37 {a) : hence the 
name “ shifting theorem 

* This is Heaviside’s unit function; the notation (21.2) will 
be found very useful in many problems. It has not appeared 
earlier in the present treatment, since problems have been stated 
hitherto in the form of, say, “ unit voltage applied to a circuit 
at f = o ”, instead of “ voltage H(i) appHed to a circuit ” : it 
must be noticed that these statements are identical only for 
t > o ; for t < o the first does not specify the voltage but the 
second makes it zero. 


82 
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To prove (21. i) we have only to notice that 

fOO fOO roO 

I — d)dt^\ e'^^^y{t — e'~^^e'^^^y(t)dt 

Jo J a Jo 

Using this theorem we can build up the transforms of 
a large number of useful step and broken functions. 

E 


0 <X t->- 

Fig. 27 (b). 

Example i. y{t) = E, o < ^ < <2^ 

= 0, t> a y 

i.e. y(t) == E{i - H(^ — ^)}, cf. Fig. 27 (b), 

j; = 5 (i — e-ap). . . (21.3) 


E 


Fig, 27 (c). 




Example 2. A half sine wave, Fig. 27 (c), 

3; = E sin o <t < 7r/a)\ 

= 0, i > tt/co /’ 

e. = E sin + E sin oj(t — 7 rjo))}l{t — Tr/oj), 
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Example 3. A “ square ” wave, Fig. 27 
3; = E, rT < f < (r + i)T \ 
= - E, (r + i)T < f < (r + a)Tj 
i.e. y = E{H(i) - 2H(t - T) + 2H(t - 2 



E I — e~^'^ 


= I tanh i-^T. 



Example 4. A '' saw-tooth wave, Fig. 
y ^k{t — rT), rT < 3 f < (r + i)T, 3 

= + + - • •) 

= A_ 

p{i ~ e“P^)' 
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Example 5. A function /(i) of period T, so that 
f{t + rT) = f{t), if r is any integer. 

_ fOO fX f 2 X 

fiP) = J ^ B-fi{t)dt = J ^e-fif{t)dt + J ^ e-P^f{t)dt H- . . . 

= ^e-v%t)dt-]re-M^e-Pif{t)dt-\- . . . 

Jo Jo 

= (i + + e-ai-T + . . .)^^e-vtf^t)dt 

= I Z^ -g T • . (21.7) 

(21-5) and (21.6) can of course be deduced from this result. 


Example 6. The transforms of the functions in Examples 
3 to 5 have all involved hyperbolic functions. A technique 
which is of great importance both here and in connection 
with partial differential equations, is that for finding the 
function which has a given transform of this type. To 
do this we expand the transform by the Binomial Theorem 
in a series of negative exponentials : for example if 

E E I — 

m = ptanh pT-=p. 

= ^(i _ e- 2 J’T)(i _e-23>T 4 . e-ipT _ . . .) 



— 3e-2i!T _|_ — . 




Then, by Theorem VI, 

y(t) = Et — 2E{t — 2T)H(^ — 2T) 

+ 2E(^ - 4T)U(t - 4T) - . . (21.8) 

the function y is shown in Fig. 27 (/) (p. 86). 


Example 7. A battery of voltage E is applied at ^ = o 
to an L, R, C circuit with zero initial conditions, and at 
t = T the battery is short-circuited. 
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Here the voltage in the circuit is E — • EH(i — T), and 
by Example i its transform is 

Then, as in (9.2), the transform of the current is 
« E(i - e'P^) 

L[(p + 

and, by Theorem VI, 

I = sin nt — T) sin n(t — T)}, (21.9) 

for the case > o. 


O 

Fig. 27 (/). 

The solution for applied voltages such as those of 
Examples 3 and 4 is found in the same way (cf. Example 5 
at the end of this Chapter). An alternative form for the 
solution in these cases can be obtained by using the ex- 
pansion formula of § 33. 

Example 8. To find the deflection y at any point of a 
light uniform beam of length /, which is clamped hori- 
zontally at the same level at its ends, a; = o and aj = /, 
and carries a concentrated load W at a? = ^. 

Using the 8-function defined in § 19, a concentrated 
load W at X = a may be treated as a distributed load 
WS(a: — a). Then, proceeding as in § 6, we have to solve 
the differential equation 

EID^y = W 8 (x - a), . 

with y = Dy = o, when x = 0 and x = L 



(21,10) 
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The subsidiary equation is, by (19.7) and Theorera VI, 
and thus, from (21. ii) and Theorem VI, 

y _ _ ^H(« - a) +-^ + ^gj., 

where 3^2 unknown constants to be found from 

the conditions y = Dy ==i o, ot x = L Solving for these 
we get finally, 

W(x-ay^, . 

y a) 

, ^ax\l-aY Wxm - a)Hl + 

2 mP 6PEI 

22. Heaviside's series expansion 

This is a useful method for determining the behaviour 
of the solution for small values of the time — the solutions 
obtained previously are equally suitable for aU values or 
the time : those of this section are useful only for small 
values of the time, but are much more quickly obtamed. 
than the general solutions. 

Most of the transforms in Chapters I and II have oeen 
found to take the form 

■,-M 


y = 


g{py 


(32. l) 


where g{t>) is a polynomial in ^ of degree w, and. fip) ^ 
polynomial of lower degree. If we divide numerator and 
denominator of (22.1) by and expand in ascendmg 
powers of ijp by the Binomial Theorem, we find an ex- 
pression of the form 


( 22 - 3 ) 


' SiP) P P 

and thus y = + ^3^ + 

by (i- 5 )> y obtained as a series of ascending powers 
of t. 
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Example i. In § 9, Example i, we had 

1== 5 

L[f.^ + (R/L)i> + (i/LC)] 

E 

L^2[i+(R/Li.) + (i/LQ>^)] 

A/ ^ ^ 


Thus 


U) 


I = 

L zU 


EfV I 

6L\LC 


Example a. , In the problem of § 14, (14.1) gives 

I_ i+(ilL,C,p^) • 

Rp{i +[(C+C,)/RCCiP] +[i /LiC^^] +[i IRCUCiP^} 


Rp\ 
Therefore, 

If we let p 


RCQiP 


+ 






I C + Q 
R R^CCi ‘ ‘ 
00 in, (aa.a), we find 

to MW 


./is -limij’W, 

3>->cosvP; p-^co 


and since, by (22*3), is the value of y when 0, we 
have proved, /or the form (22.1) ofy^ the following Theorem. 


Theorem VIL If y{p) is the transform of y{t)^ then 

lim 3;(^) = lim ^37 (jP). . . . (224) 

p-T^co 

This result is in fact true under fairly general conditions. 
The following theorem, of the same type as (224), may 
often be used to determine the behaviour of a function for 
large values of the time : 


Theorem VIII. If y[p) = f{p)jg{p)y where f{p) and 
g{p) are polynomials in p, the degree of f(J>) being less than 
that of g{p)y and if the roots of g{p) = o are either zero or 
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have negative real parts then if y{t) is the function whose 
transform is y[p) 

\im y{t) == \\m py{p) , . (22.5) 

>00 5 )— >-0 

and^ under the same conditions^ 

poo 

y{t)dt == lim y{p) . . (22.6) 

J 0 

The result follows immediately from the formulae of 
§ 2 for expressing y{p) in partial fractions. An application 
of (22.6) has been given in (11.8). 

23. Network theorems 

Many of the general theorems of steady state network 
theory have immediate generalizations in terms of Laplace 
transforms. 

Here we mention only Thdvenin’s theorem because of 
its importance in simplifying the algebra in complicated 
networks. A special case of this theorem, which is ade- 
quate for many purposes, may be stated as follows : 

Suppose AB and CD are two pairs of terminals of a 
network ; let Z be the generalized impedance of the net- 
work looking in at CD with AB short-circuited ; also let 
V be the transform of the open -circuit voltage drop across 
CD due to voltage V(^) applied for ^ > o at AB with zero 
initial conditions. Then the transform of the current I 
in an impedance connected across CD, due to voltage 
V(^) applied at AB for > o with zero initial conditions, 


is 


1 = 


Z + 'S’o 


(23 -i) 


As an example, suppose the network is a filter circuit 
of 7n, T sections [Fig. 18 (a)\. In this case the quantity 
Z defined above is, by (16.10) with r = o, 

Tj — z sinh 0 tanh mQ . . (23.2) 


’••‘This Theorem is frequently stated in the literature with 
no conditions on yip'). Clearly it does not hold, for example, 
if yif) — sin tat. 
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and the quantity v defined above is, by (i6.n) with 


y = 7/r — I 


= V sech md 


Thus, by (^3.1), the current I in an output impedance 
Zq is given by 

i = _ . 

z sinh 0 sinh mO + ;8:o cosh md 

This agrees with (16-7), but the calculation is simpler, 
since Th6venin’s theorem allows the result for any terminal 
impedance to be deduced from the solutions for open and 
short-circuited terminations and these are often known 
and in any case are e^y to calculate ab initio, 

24. The superposition * theorem 

Theorem IX. If y-iip) is the transform of yft)^ 
yziP) transform of yft)^ then yx{p)yfip) is the trans- 

form of 

\ yl{^')y^{t — t')dt' = [ y^{t — t’)y^{t’)dt' . (24.1) 

Jo Jo 

The proof f is not particularly difficult, but is omitted 
here as it depends on the transformation of double integrals, 
and involves only routine pure mathematics. 

* Also referred to as the Composition, Convolution, or Faltung 
Theorem, and sometimes as BorePs or as DuhamePs Theorem. 

t C. and J., §33 ; Doetsch, loc. cit,, p. i6i. The essential 
steps, which require justification and some conditions on the 
functions, are as follows : 

( C30 *00 

where the double integral is taken over the quadrant u > o, 
V > o. In this put u v — v = t', then t' can take all values 
up to t, and t all values up to co , and the double integral becomes 


Jo Jo 
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As a first application consider the differential equation 
(D^ + a^)y = f(t), t>o, 

to be solved with y = D3; = o, when t = o. The sub- 
sidiary equation gives 

m . ■ 

Thus by Theorem IX 

y = ^ f f{t') sin a{t — 

0 

Clearly, any system of ordinary linear differential 
equations with constant coefficients and with arbitrary 
functions on the right hand side can be treated in the same 
way. ^ ^ ^ 

From the present point of view, the chief use or 
Theorem IX is the following : suppose we require the 
current I at some place in a network due to voltage V(^) 
applied at some point with zero initial conditions.* In "the 
usual way we find 

1 = 

where z{^') is a known function of p. Then if x(t) is the 
function, found in the usual way, whose transform is 
I Theorem IX gives 

I = I - t')dt' = I* V(f - t')x{t')dt' (24.3) 

Thus for any applied voltage V, which may be so com- 
plicated that we cannot evaluate its transform, or which 
may be given numerically, for example by an oscillogram, 
we can always express the current by a definite integral 
which may be evaluated numerically. In the same way, 
if the current in any part of a circuit is known, the current 
in any other part can be determined. 

* For non-zero initial conditions there will be ad^tional terrns 
on the right of (24.2) which can be dealt with as in Chapter I . 


Y(p) 

z{py 


(^ 4 - 2 ) 
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For example, for voltage V applied to an L, R, C circuit 
with zero initial conditions, we have as in (9.2) 

T 

L[(p + 

and thus, if > o, 

I = cos nV — yb sin nt'}e~‘^^* df (24.4) 

The function x{t) used above whose transform was 
ilz{p) is, by (24.2), the current due to a unit impulsive 
appHed voltage, V = S(^), = i. This is an interesting 

transient in its own right, and the present result shows 
that if it is known, the current due to any applied voltage, 
V, may be evaluated as a definite integral. The other 
simple and interesting transient in the circuit is the 
current due to unit applied voltage, V = i, which, by 
(24.2) has the transform 



If I3L has been found as a function of t from (24.5), we 
may express the current due to any applied voltage V as 
an integral involving Ii, by writing (24.2) in the form 

M 

Now, by Theorem IX, the function whose transform is 

V(p)/l>«(/>)] is 

J\(i')V(i - . . ( 24 - 7 ) 

Then it follows from (24.6) and Theorem III that 

since (24,7) vanishes when t — o. Clearly various trans- 
formations of (24.8) are possible. 

Finally, it should be remarked that the importance of 
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Theorem IX is in the study of difficult or arbitrary func- 
tions, or functions given numerically by oscillograms ; if 
V is a simple exponential or trigonometric function, it is 
simpler to evaluate I directly by the methods used earlier 
than to evaluate the integrals (2^4.3) or (24,8). 

35. The inversion theorem 

This is a general theorem* which gives y{t) from 
subject, of course, to conditions on y{t) orj;(^). Omitting 
these conditions, it states that if 

yi^) = ^ R(/>)>o, . (25-1) 

j ry+ico 

then y{t) == — .1 e^*y(AyA, y > 0. . (25*2) 

27rtJ y-ico 

In (25.2), A is a complex variable,' and the integral is to 
be interpreted and evaluated by the theory of functions 
of a complex variable. If this theory, and in particular 
the calculus of residues, is known, (25.2) provides an 
important method for finding the function which has a 
given transform. For the problems of Chapters I and II 
which led to ordinary linear differential equations, this 
gives an alternative method for finding y{t\ which, how- 
ever, is no shorter than those used earlier, so that for such 
problems there is no necessity to learn the technique of 
complex variable. For a thorough study of the applications 
of the method to partial differential equations, the use of 
the inversion theorem is necessary. 

From the point of view of the historical development of 
the subject the correspondence between (25.1) and (25.2) 
is of great importance. Carson *[• regarded (25.1) as an 
integral equation from which the solution of the problem 
could be obtained. Bromwich, and to a certain extent 

* C. and J., §§ 29, 30. Doetsch, Zoc. cit,, Kap. 6. The result 
is sometimes called the Fourier-Mellin Theorem. It is related 
to Fourier’s and MelHn’s inversion formulae, cf. Doetsch, loc. ci£. 

t Electric Circuit Theory and Operational Calculus (McGraw 
mi, 1926). 
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Wagner, Jeffreys, and McLachlan, regard an integral of 
the type (25.2), with j?(A) found from the usual subsidiary 
equation, as the solution of the problem. The Laplace 
transformation, and its inversion theorem provides the 
link between these two ideas ; this point of view is largely 
due to Doetsch. 


26. The connection with Fouriefs integral theorem and 
ParsevaVs theorem 

Fourier’s integral theorem* states that, subject to con- 
ditions on y{t\ if oj is real and 

Y(a)) = r . . . (26.1) 

J —00 

then e'~^^^Y(a))dcx). . (26.2) 

—CO 

The simplest of the essential conditions on y(t) is that 



be convergent. 


(26.3) 


Parseval’s theoremj* states that, roughly provided these 
integrals exist, 


f “ I y(t) 1 2 if = — f " I Y(a)) 1 2 io). . (26.4) 

J —00 27r J _oo 

For comparison with the problems previously discussed, 
suppose that y{t) is real, and is zero J when t <0, Then 
(26.1) and (26.2) become : 


♦ For a proof see Carslaw, Fourier Series and Integrals, edn. 3, 
§ 119 ; the result (26.3), (26.1) follows from his on using the 
exponential form for the cosine. See also Titchmatsh, Theory 
of Fourier Integrals y §§ i.i, 1.2. 

tTitchmarsh, loc. cit.y §2.1. The theorem is often referred 
to as Rayleigh’s theorem, and its use in physical applications is 
due to him [PAiZ. Mag.y 27 (1889), 466]. 

J It is frequently stated that the Laplace transformation method 
contemplates only functions of this type. In the method as 
developed in Chapters I and II the problem was to find the solu- 
tion of a differential equation for f > o, which satisfied certain 
conditions at t = o, and in fact the solutions found would also hold 



THEOREMS ANP THEIR APPLICATIONS 


if Y{co) = • ( 2 ^- 5 ) 


then * 

-R(”e-*-“’‘Y(a))<ia) = y(t), H t > o 
Jo 


(26.6) 

and 

-Rf“e-‘'‘“*Y(a))^f£o = 0, if « < 0 

■TT Jo 



provided 

rco 

\ 1 y(f) 1 dt is convergent. 

J 0 

• 

(26.7) 

To deduce (26.6) from (26.2) in this case, we 
that the right-hand side of the latter is equal to 

notice 


-i(V<‘“‘Y(a.)ia) + i-("e^"‘Y(— cu)<ia>, 

27 rJo 277 Jo 


and, since is real, Y(— co) is the conjugate of Y(cu), and 


the result (^6,6) follows. 

Also, in this case (26.4) becomes 

ri:);(i)12rfi = iriY(a>)l® ifco. . (26.8) 

JO ^Jo 

Y(a)) in (26.5) is called the Fourier transform of y{t). 
Since (26.5) is precisely the same as the equation defining 
the Laplace transform y{p) of y(t)f with p replaced by 
(— io)\ it follows that, ifY{<jS) exists, 

Y[aj) — y{— i(jS)y . . ( 26 . 9 ) 

Thus the Fourier transform is the degenerate form of 
the Laplace transform if the latter exists when p is pure 
imaginary, also (26.6) is the degenerate form of (25.2) if 

for t<o; for example, the solution of (D 1)3^ — o with 
y = I when t = o is found to be y — ^ e’-^ for either i > o or 
t < o, though we are only interested in the fornner caac. In 
this section and in § 25 there is a difference in principle (though 
not in the results obtained for t > o) since a solution ia sought 
which is zero for t < o \ for example the solution of the differ- 
ential equation above, found from (25.2), would be :v(0 o, 
t<o] y{t) = i, i =: o ; y{t) = e"*, t > o. 

* R is written for ** the real part of ”, 
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y can be taken zero. The Laplace transform exists for 
many functions for which the Fourier transform does not 
(for example, i, or sin t) but, if the latter does exist, it 
can be found from (26.9) and previous results. 


For example,* if y{t) = e“"%in 

iS 


so that 

Y(co) = j 7 ( — tco>) 


y{p) = 


iP + a)2 + 


(a-ico)2-f;82 (a2+;82. 

In this case (26.8) gives 

rco jS2foo / 

Jo 7rJo(a® + P“ — < 


dcxj 


■ + 4a 


. (26.10) 


- a>^) — zicf.(x> 

. (26.11) 


From the point of view of electric circuit theory the 
great importance of Fourier’s integral theorem lies in its 
physical interpretation. The Laplace transform has no 
obvious physical significance ; it is a purely mathematical 
device for finding a solution — this is not altogether a 
disadvantage as it implies that any difficulties can be form- 
ulated and treated mathematically, while it was on points 
of interpretation that the old operational method was 
obscure. On the other hand (26.6) expresses y{t) as a 
combination (that is the integral) of harmonic terms 
each with its Own amplitude Y{o>)do}^ in the same way 
that a periodic quantity is expressed by Fourier’s theorem 
as a sum of fundamental and harmonics, each with its 
own amplitude. 

Now suppose a voltage y{t) is applied to a circuit with 
zero initial conditions, and the current I in any branch of 
it is calculated, then, proceeding as in Chapter II we find 

• • • (^^-^ 3 ) 

* Fourier transforms of many functions of this type, with 
numerical illustrations and applications, are given in Burch and 
Bloemsma, Phil. Mag. (6), 49 (1935), 480. 


f 
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where z{p) is known and depends on the circuit. ‘Then 
if y{t) has a Fourier transform, this will be j 7 (— iw), by 
(26.9), and (26.13) gives for the Fourier transform of I(^) 

= 

Thus, by (26*6) 

This expression has a simple physical interpretation, 
since iw) is the steady state current in the circuit 

due to voltage Thus (26.15) may be interpreted 

by the statement that the current due to voltage y{t) is 
obtained by superposing the effects of all the component 
frequencies of 3;(t), each with its appropriate amplitude 
y(— ioi)dai. That is, we have an extension of the results 
of steady state alternating current theory to non-periodic 
applied voltages, although, for the solution of specific 
problems, this is inferior to the Laplace transformation 
methods described earlier. 

One application of (26.14), however, is new and of 
great importance. Applying Parseval’s theorem (26.8) to 
(26.14) gives 

^ dw . . (26.16) 

Jo ttJo l5'(-- 2 a>) 

The quantity on the left gives the energy absorbed in 
the branch of the circuit, alternatively it is the total re- 
sponse of a square law detector, which is of considerable 
practical interest, and (26*16) gives this as a real infinite 
integral without having actually to evaluate the current. 

For example, if voltage is applied at f = o to an 
inductive resistance R, L we have 

I dco 

ttJo I (R — L2a>)(oc — zc/>) 1^ 

I doi 

“^Jo L^(a>^ “b ^^)(a>^ -f- (X^) [continued overleaf 




,58 


THE LAPLACE TRANSFORMATION 


I ri ^ 1 CD I _ £D“1'^ 

I 

2LW(a+W’ 


where JS = R/L, 

In more complicated problems the integral (26.16) may 
be evaluated numerically, and approximations can be 
introduced, for example, with a selective circuit, [^(/cd)! is 
very large except in a certain range of values of cd, and 
(26.16) need be evaluated only over this range. 


examples on chapter III 


I. If y(t) is I, 2, 3,.. .. in the regions o < f < T, T < f < 2T, 
2T < i < 3T, and so on (a staircase ” function) 




' p(i — 

2. If y(t) =1, o < f < fi ; 3;(f) = 0, ti<t<T; and this 
is repeated with period T (a repeated pulse, or the difference of 
two staircase functions) 

. I —e~Ph 
“ ^>(1 - 

3. If y(t) — I sin out |, (full wave rectified alternating voltage) 
6 )(i + 


y(P) ■■ 


+ a)^)(i — 


4. A uniform beam is clamped horizontally at both ends. It 
carries a distributed load tv (constant) per unit length in o < x <il 
and zero in il < x < L Show that the deflection at any point is 

I3zvlx^ iiwPx^ tvx* w(x — iOVir 

192EI 384EI 24EI 24EI 

5. If the square-wave voltage of § 21, Ex. 3, is applied at . 
t — o to an L, R, C circuit with zero initial current and charge, 
show that the current in the circuit at time t = zsT + i', where 

s is an integer and t' < T, is 
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using the notation (9.3) and assuming > o. Deduce that 
in the steady state due to a square wave voltage, the current at 
time after the beginning of a cycle is 

2E sin ni' — sin n(T — t') 

nLi I + cos wT + 

6. In the problem of § n, Ex. 2, show that for small values 
of the time 

T — ^ -L \ 

^ L V 6 60 ‘ 


7. In the circuit of Fig. 7, an unknown voltage is applied at 
S at t = o with zero initial currents and charges. If the voltage 
drop across the first condenser is show that the voltage drop 
across the second condenser is 


n f vit*) sin w(f — t')dt\ 
J 0 


8. If y{t) = t > o, show that its Fourier trans^ 

form y(<w) is 


(a — — io))' 


Show that the total energy absorbed by an inductive resistance 
D, R from this voltage is 

R(j 3 - ay( 0 L d- jg + /^) 

2L^ajS/A(a + p)(fi H- -h a)’ 


where /x = R/L. 



CHAPTER IV 


PARTIAL DIFFERENTIAL EQUATIONS 
27. Introductory 

In the preceding chapters it has been found that problems 
leading to ordinary linear differential equations could be 
solved with only the elementary apparatus of §§ i-3- 
the study of partial differential equations it will appear that 
it is still easy to find the Laplace transform of the solution, 
but, since this usually involves more complicated functions 
of p than those considered in Chapter I, it is rather more 
difficult to derive the solution from it. The Table of 
Transforms can of course be enlarged, but the most power- 
ful general method is the application of the inversion 
theorem, § 25, combined with the use of contour integration 
to evaluate the complex integral which occurs in (25*2). 
For a complete study of the subject, particularly in its 
theoretical aspects, a knowledge of this technique is essential ; 
but it is possible to go far without it, for example most 
results on finite transmission lines, and some on infinite 
lines, can be obtained, in some cases rigorously, and in 
others with the proviso that for a complete treatment the 
complex variable would have to be used, although the 
actual process of writing down the solution would be 
precisely the same as that used here. 

The differential equations for a uniform electric trans- 
mission line will be found in § :a8, and, for shortness, only 
these will be studied, although many of the problems 
considered are problems in wave motion or conduction of 
heat in one dimension, with a slightly altered notation. 

28, The dijferential equations of the uniform transmission 
line 

We suppose the line to have resistance R, inductance L, 
capacity C, and leakage conductance G, per unit length. 

ICX> 
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Also let I be the current in the direction of increasing x 
at the point x of the line, and let V be the voltage 
drop across the line at this point. I and V will be func- 
tions of both distance x along the line and time t ; to 
emphasize this they will sometimes be written i) and 
V(^, t). 

To find the differential equations satisfied by the current 
and voltage in the line, consider the points A at x, and 
B at The current and voltage at A will be I and V, 

and at B they will be I + and V + ^Sx, respectively. 

The portion of the line between A and B can be represented 


LSx R8x 




i 


Fig. aS. 






c + Sx 


by the equivalent circuit of Fig. 28. Here the voltage 
drop from A to B over the series portion of the line is 

— and the current in it is I + ; neglecting 

Tix ox 

terms in (8a:) ^ this gives 


^^ + R1)Bx=-^Sx 


It 


(zS.i) 


Also the voltage drop over the parallel portion HK of 
the line is V, and the current in it is — ^So?, thus 


( -I- GV^Sx == - ^S;c. 
\ It J 


()L 


. (28.2) 
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From '(28.1) and (28'.2) it follows that the differential 
equations for I and V are 

’=S+°''=-S ■ ■ 

These are a pair of simultaneous partial differential 
equations for I and V, We wish to solve them for ^ > o 
with given initial values of I and V, which we shall write 
!(:«) and V(«) to emphasize that they are functions of x. 
Explicitly these initial conditions are 

I(a;, t) ^(^^[;), 2 L 3 t^o for fixed . (28.5) 

V(a;, t) V(^), as ^ o for fixed x. . (28.6) 

There are also boundary conditions to be satisfied at 
the ends of the line; these depend on the particular 
problem. For example, if the line is short-circuited at 
X = we must have V = 0 for x = a and all i > o ; 
if there is open circuit at a? = a, the condition is I = o for 

X == a and all ^ > 0 ; if the voltage is prescribed at a: = a, 

V(a, t) must be this prescribed function of i for f > o ; 
if the line is joined to some other impedance at x = a, 
there will be a condition involving both the current and 
voltage at x = a, 

29. The solution of partial differential equations by the 
Laplace transformation 

We consider first the system of equations (28.3), 
(28.4) which have to be solved for ^ > o with initial con- 
ditions (28,5) and (28.6),' and with boundary conditions 
to be specified. Proceeding * precisely as in Chapter I, 
§ 4, we multiply the equations by 6"^^ and integrate with 
respect to t from o to 00, Then, just as in the proof of 
Theorem III, 

* As in Chapter I assumptions about the function I are implied, 
and, strictly, the solution should be verified. 
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i)T"“+ P re-®‘I(*, t)dt 
Jo L Jt=o Jo 


- +^i 


(29.1) 


where I is now a function of x and py and, when we wish 
to emphasize this, it is written l{Xy p). 


where we have assumed that the orders of integration, and 
of differentiation with respect to the parameter Xy may be 
interchanged. 

Using these results, and the corresponding ones for V, 
we obtain from (28.3) and (28.4) and the initial conditions, 
the “ subsidiary equations 

(Lp + K)l{x,p) = - + Lt(*) (29-3) 

(Cp + G)?(*, p)^- + CV(*) (29.4) 

Eliminating I gives the ordinary differential equation 
for . 

g_j^ = Lg-qi^ + R)V, . (»5.s) 

where = (Lp + R)(C/) + G) . . (29.6) 


When V has been found, I is given by (29.3)) that is, by 
[ . Lt{x) 


1 = - 


Lp dx Lp + R 


(29-7) 


It appears that the effect of the Laplace transformation 
with respect to t is to reduce a partial differential equation 
in A? and ^ to an ordinary differential equation in x. This 
has to be solved with boundary conditions which are the 
haplace transforms of those originally prescribed — ^they are 
hest written down as they arise in specific problems. 

The quantity q which appears in (29.6) is, in the general 
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case, the square root of a quadratic function of but in 
three important cases it takes a simpler form. These are : 

(i) The lossless line ”, R == G = o, for which 

q = pjc, where c = (LC)“^ . . (29.8) 

(ii) Heaviside’s ‘‘ distortionless line ” in which 

R G 
L - C ~ 


and thus - ^ + p)/^, where c = (LC)""^. . (29.9) 

(iii) The ideal submarine cable, L = G = o, for which 

5 = -v/(RCp)- - ■ (29- 10) 

Another partial differential equation of importance is 
the equation of linear flow of heat {v temperature, /c 
diffusivity) 


I 'bv 


(29.11) 


to be solved with v = »(*), when t = o. The subsidiary 
equation, found as above, is 

d^v p_ I. . 

. . (29.12) 


and thus is of the same type as (29.5) in the case (29.10) 
ofL=G = o. 


The equation of wave motion in one dimension is 


(29-13) 


where y is the displacement. This has to be solved with 
y = y(^)> ^ when t = o. 


The subsidiary equation is 


d^y 




4 v( ' 


*• Tr/ \ 
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which is of the same type as (29.5) for the case (29.8) of 
R = G = 0. 

When the subsidiary equation (29.5) with its boundary 
conditions has been solved,* we know the transform 
P) solution V(^, t). For problems on semi- 

infinite lines we rely on the Table of Transforms to find 
the solution from its transform. For problems on finite 
lines two methods are available : (i) an extension of the 
expansion theorem (2.8), which gives the solution in the 
form of a trigonometric series ; or (ii) an expansion in 
terms of the solutions for the semi-infinite line, which has 
a physical interpretation in terms of direct and multiply 
reflected waves. Both these methods will be used, treat- 
ing the latter first. 

30. The semi-infinite line 

We consider the problem of the semi-infinite line, 
A! > o, with zero initial current and charge. First, suppose 
the end ai = 0 to be maintained at constant voltage E for 
£ > 0. 

In this case the subsidiary equation (29.5) becomes 


= X> 0 . . 

• (30.1) 

The boundary conditions to be satisfied at the ends of the 
line are 

V = E, x = 0, t> Q, . . (30*2) 

V to be finite as « 00 

• (30-3) 

The transforms of (30.2) and (30.3) are 


II 

W 

II 

0 

• (304) 

V to be finite as a; 00. . 

• (30-5) 


We could, as in § 6, solve it by the^ Laplace transformation 
method, but, for the simple problems arising here, it is just as 
easy to write down the particular integral and complementary 
function in the usual way. 
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The general solution of (30.1) is 

V = + Be®*®, 

and (30,5) requires B = 0, and then (30.4) gives A = 

Thus, using the value (29.6) of 

^ . (30.6) 

P P 

In the general case the transform involved in (30.6) is 
quite different from anything previously encountered, and 
in fact V turns out, after a difficult calculation which will 
not be given here, to be a complicated integral involving 
Bessel functions (cf. § 35). In the important special cases 
mentioned in § ^9, however, V can be evaluated simply. 

For the lossless line, R = G — o (30.6) becomes 

^ . . . (30.7) 

P 

where c = (LC)-^. . . . {30.S) 

Then it follows from (a i.i), with y(t) = E, that 

V = EH[f-(*/c)], . . (30.9) 

that is : V is zero up to the time xjCy at which a wave 
travelling with velocity c from the origin would arrive ^ at 
the point and V has the constant value E at this point 
for greater times. 

For the distortionless line (29.9) we have 

■p 

"V* = , . , (30.10) 

P 

and ■ V = Ee-W«H(^ — xjc) . . • (30.1 1) 

In this case the disturbance is still propagated with 
velocity c, but there is attenuation as we move along the 
line. 

Some problems on the submarine cable, L = G = o, 
will be studied in § 35. 

Suppose, now, that instead of applying a constant 
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voltage E at A! = 0, we apply a sinusoidal voltage E sin 
The only modification is that in place of (30.6) we have 




Eo^ 






(30.12) 


In the case of the distortionless line, ^ = (^ + p)/^> 
this becomes 

• • (30-13) 

and V — sin m{t — (30-14) 


Still more generally, for any applied voltage, /(i), we 
have in place of (30.6) 

V=/(p)e-«*. . . (30.15) 

For the distortionless line this becomes 

‘V =^^)e“-(3>+p)®/c, , (30.16) 

and, again using (21. i), 

V = . (30.17) 

This may be interpreted by the statement that the 
voltage at x is zero up to time xjc, and subsequently follows 
that at jc = o with a time lag of xjc^ and a reduction in 
magnitude by the factor 

In the next section the application to problems on finite 
lines of the solutions found above will be studied. 

31. The travelling wave solutions for the finite transmission 
line 

As an example of the method we consider the problem 
of the finite transmission line 0 < ix; < Z, with zero initial 
current and charge. The end ac = 0 is to be earthed, and 
the end a? == Z maintained at constant voltage E for Z > o. 
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The subsidiary equation (29*5) becomes in this case 

^-q^Y = o, o<x<l. ' . ( 31 - 1 ) 

The boundary conditions are 

V = E, « = Z, f>o, . . (31-2) 

V = 0, X = o, o. . . 


The transforms of (31.2) and (31.3) are 

Y^Elp, x=^l, . ■ ( 314 ) 

Y = o, x = o. . . (31-5) 

The solution of (31. i) which satisfies (31.4) ( 3 ^-S) 

= . . (3X.6) 

p smh ql 


Here q is defined by (29.6) so that 'y* is a complicated 
function of p. One method of proceeding is suggested 
by the device of expanding in a series of negative ex- 
ponentials which was used in § 21, Example 6. We write 
(31.6) in a form involving negative exponentials, and 
expand the denominator by the binomial theorem. This 
gives 

'~p * 
p' 

— g-gf(J+a5)_j_^-fl(3Z-aJ) — ^-a(3Z+a)q_ ^ ^ (3 1 - 7 ) 

In the special case of the lossless line, for which, by 
(29.8), q = pjc^ (31*7) becomes 

tjr = ^^q—t(1—x)Io — g— jJ(Z+a;)/c g~p(3Z~a;)/o — ^ ^ 



PARTIAL DIFFERENTIAL EQUATIONS 109 
And, using (21. i) 

• (3 1 -9) 

The graph of the function (31.9) is shown in Fig. 29. 
The voltage at x is zero up to time {I — x)lc, at which a 
wave travelling direct from the end x = I would reach the 
point X, The voltage then has the constant value E up 
to the time (/ + x)/Cy at which a wave travelling from the 
end = / and reflected back from the end x = o would 
arrive. From this time up to the time of arrival of a 
twice reflected wave, it has the value zero, and so on. 

Direct Once Twice Thrice reflected 







E - 



1 t3t+x)/c 


Fig. 29. 


Returning to the general expression (31.7), it is seen 
that the terms of this are all of the form which 

appeared in (30.6) in the corresponding problem for the 
semi-infinite line. Thus, if the solution for this case is 
known, the solution for a finite line can be expressed in 
terms of it. 

As a second example, we consider a lossless line of 
length Z, initially charged to constant voltage E. At 
i = o the end x = o is earthed, the end x = I being left 
insulated. 

In this case the subsidiary equation (29.5) becomes, 
writing c = (LC)'"^, 

dx^ c"’ 


• (31*10) 



no 
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which has to be solved with 

7 = O, X = o, . . (31-11) 

1 = 0, * = /. . . (31.12) 

Also it follows from (29.7) that 

The solution of (3 1 , 10) which satisfies (3 1 . 1 1 ) and (3 1 . 12) is 

V==--® cosh j.(l - x)lc , s 

p p' cosh pljc • ■ • ‘ 

]g ^-33a/cJj; ^ 

'~'pV I -I- e-ipiio / 

E 

= — .Jl — ^-pxlo — ^-p(2l-x)lo ^g-3){^l+oa)lc , .J. 

V 


Thus J>.-h(.-”)-h(<-V") 

)+.. 


. +“(‘ 

Also by (31.13) 


2I + X' 


(31-15) 

(31-16) 


C'^ si nh x)lc 
cosh pljc 

The most interesting quantity is the current at = o. 
At this point 

• _E //C\ i 
p \ \L/ 1 + 

(31-17) 

that is, a square wave of amplitude E-\/(C/L) and period 
(4I/C). 



PARTIAL DIFFERENTIAL EQUATIONS 


III 


32. The finite line with terminal impedances 

Suppose that at = Z, where the voltage is Vj and the 
current is Ij, the line is terminated by impedance Zi^ 
Then, for zero initial conditions in this impedance, the 
boundary condition 2X x — I 2c=o 
is just the subsidiary equation y Vg 

for the impedance Zi^ namely 

arjli = ■?,. . . (32.1) ZoIJ ■ 

Similarly, if voltage v is yie 

applied through impedance X ^ 

5^0 to the line sit x = o ^ 
writing Vg for the voltage at 

the point ix; = o of the line and Ig for the current there, 
the boundary condition at x ^ o, for zero initial con- 
ditions in the impedance Zq is 

Zo^o = w — Vo. . . . (32.2) ■ 

As an example, suppose that a lossless line of length Z 
is initially charged to voltage E, and that at ^ = o it is 
discharged through a resistance Rg at ^ = o, the end 
X =: I being insulated. In this case th^subsidiary equation 
(29.5) becomes 

0<*</, . (32.3) 

where c ~ (LC)"^'. At == Z, we have 

1 = o, = /, . . . (32.4) 

and at = o the boundary condition (3^.^) is 

Rfllo = - Vo- • • • (32-5) 

Also we know from (29.7) that 


T- 

hpdx' ■■ 


and thus (32.5) becomes 
Rq^V 

Lipdx 


-j- "V == o, X = o. 


(32.7) 
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The general solution of (32.3) which satisfies (32.4) is 

^ = 5 + A cosh p{l ““ x)jc. . (32^8) 

Substituting this in (32:7) gives 

sinh — + A cosh ^ + 5 = 0. 

Lc c * c p 

Using this value of A in (32.8), we find 

^_Er Lc cosh p{l^x)/c \ . . 

p\ Kq smh {pi I c) he cosh, {pi I c)) 

The current is obtained from (32.6) which gives 

T = _ Esmhj>(/-jg)/c / -v 

P{Rq sinh {pl/c) + cosh {pHc)Y 

where = Lc = ^(L/C). For shortness we evaluate 
only the current at = o. In this case, expressing (32.10) 
in negative exponentials, and writing a=(Ri— Ro)/Fi+I^o)> 
we find 

T - E(i 


}• ( 32 - 9 ) 


1 = - 


^(Ri + Ro)(i + «e-23>*/o) 


= - W+rA ' + <“ + 


Therefore 


' ^ I E(« + i\ .( 

Ri + Ro^Ri + Ro"V 


(32- 1 0 


+ (32.12) 


so that the current has a series of constant values, diminish- 
ing in magnitude. If R^ = R^ = ^(L/C), (the character- 
istic impedance), there is a single pulse of current — E/2R0 
for time zljcy and zero current subsequently. 
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Clearly any terminal impedances can be discussed in 
the same way, but since the coefficients of the later ex- 
ponentials are usually rather complicated functions of p, 
the method is most suitable for use when the solution is 
needed for fairly small values of the time. 


33. The expandon formula for finite Urns 

In chapters I and II the transforms of the solutions of 
problems on ordinary differential equations were usually 
found to be of the form 







gipy 


( 33 -:^) 


where f{p) and g{p) were polynomials in p. The solution 
was found from ^(p) by expressing (33.1) m partial 
fractions, or by the equivalent technique of Theorem II, 
for example, by the result (2.8) that if are the 

zeros of^(^), and these are all different, then 

• • (33.2) 


In the problems on finite lines in this chapter, rather 
more complicated transforms have appeared. The new 
feature of these is the occurrence of quotients of hyper- 
bolic functions of q» For example, the expression 


_ cosh j)(/ - x)lc 
p cosh. pi jc 


(33-3) 


appeared in (3 1.14). Since cosh z and sinh z can be 
represented by the infinite products * 


sinh s _ ^ 

* Carslaw, Plane Trigonometry ^ edn. 3, § i68. 


(334) 

(33-S) 


8 
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a quotient of hyperbolic functions such as (33.3) may still 
be regarded as being of the type 

■ ■ ■ < 53 . 6 > 

except that nov^ f{p) and g{p) may have an infinite number 
of factors, instead of a finite number as in (33.1). Then 
if ^2, . . . are the zeros of g{p), and if these are all 
different, it is plausible to assume that (33.2) can be ex- 
tended to this case : that is, if ^[p) is given by (33.6), 

Just as in § 2 (33.7) apj)lies only to the case in which 
g{p) has no repeated zeros ; if it has, the corresponding 
extensions (2.19) or (2.20) must be used. 

The result (33.7) will be referred to as the expansion 
formula ; it has not been proved, and is merely suggested 
as a possible extension to the case of an infinite number of 
zeros, of the result (2.8) which was proved for the case 
of a finite number of zeros. * That is precisely the approach 
which was, used by Heaviside. The justification for the 
use of (33.7) is as follows : the rigorous pure mathematical 
method of finding from (33.6) consists of the applica- 
tion of the inversion theorem (25.2), followed by contour 
integration and a certain amount of pure mathematical 
discussion, and, in the case of finite lines the final step in 
this procedure is equivalent to the use of (33.7)- Thus it 
may be taken that there is a solid theoretical foundation 
for results on finite lines (or on wave motion, or conduction 
of heat, in a finite region) obtained by the use of (33*7) J 
this remark does not apply to semi-infinite lines for example 
(33.7) must not be applied to (30.6). 

* Any applied voltages must be combinations of the functions 
of Table I, and the lines must not be terminated by their char- 
acteristic impedance. 

t Or to lines terminated by their characteristic impedance — 
in .this case the transforms encountered are of the same type 
as those for semi-infinite lines. 
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As a first example of the use of (33.7), we apply it to 
find (j){t) from (33.3). First the zeros of the denominator 
must be found. This could be done by putting* z = pljc 
in (33.4), but, from the point of view of the subsequent 
calculation, the following alternative method is better : 
clearly == o is a zero, and the others are 


= i 


{zn + i)i'rrc 


w == 0, I, 2, 


•. (33:8) 


since, for these values of^, 

coshL^ = = cos (2^ 

C y 2I ci 2 

To apply (33.7) to (33.3) we need 

^(p cosh — = cosh — — sinh — . 


Also, since cosh {pljc) = o when p has one of the values 
(33-8)) we find from (33.9) 




^ {zn rl- i)iTr (aw + 


:( 2 n+l)wro/ 2 Z 


= K— ( 33 -II) 

Using these results, (33.7) gives from (33.3) 

^(i)=:+ f 51^1 coshr^2£±i[)i?^_«) -|,(2n+iW 

L z \ l/J -1 

+ Conjugate! 




(zn+l)irct 


^ I . (2«+l)7rUC (zn-\- l)TTCt 

A 7 i — T sin cos P — , 

^(zn-\-i) 2I zl 


( 33 -ia) 


* (33*4) and (33.5) can be used to show the actual factors of 
numerator and denominator : this is often useful in cases where 
a factor can be cancelled between them, or where some exceptional 
factor arises. 
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Il6 

where the term i has arisen from the zero p — o oi the 
denominator. 

Using this result in (31.14) gives 


__ 4E^ I . {zn-Vi)'nx (2W + 

V — — > } ; — V sin — cos i — . 

TT “ {zn + i) 


(33-13) 


A solution of the same problem, expressed in terms of 
successively reflected waves, was given in (31.15). The 
connection between this and (33.13) is of course that at 
any point a;, the solution (31.15) is a periodic function of 
t with period 4//^, and if this function is expressed as a 
Fourier series, the result is (33.13). 

As another example of the application of (33.7) we deter- 
mine V for the first problem of § 31 from its transform 
(31.6), namely 

^ Esinhg^ . 

^.sinhgr 


considering various forms of which was defined in 
(29.6) by 

= (Li. + R)(C^ + G). . (33.15) 

The first point to notice is, that if we use the infinite 
product (33.5) for sinh z in (33.14), we get 


^ ^ E qx[i + (gV/77^)][i + {q^x^lz^ 7 T^)][ , , . 

P • ql[i + {qH^ni + {qn^lz^^^)][ . . . 

<i+(^>^)(Lj) + R)(Cj> + G)][ . . . . 

P * /[I + {l^l^^)(Lp + R)(C^ + G)][ . . . 

It follows that although V, as given in (33.14), appears 
to be a function of q (which is the square root of a function 
of ^), it in fact involves only q^, and thus is a function of p 
which does not involve square roots. The same result 
holds for all transforms which arise in problems on finite 
lines. 

The generd case of (33.14) will be studied in Example 3 
below ; special cases are given in Examples i and 2. 
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Example i. The ‘‘lossless” line, q ^ pjc , where 
c = (LC)-i 

In this case (33.14) becomes 

,-5. Esinh*;w/^: , v 

One zero of the denominator is ^ = 0, and the others are 
p=:±tn 7 rc/l, w=i, 2'... . (33.18) 

. , fimrc h . . 

since sinh -! "~j- • i sm n-n* = o. 


z sm n-TT = o. 


Also r -^(p sinh — cosh — 

\_Clp\ ^ / Jp^irmoll ^ Jp^imcjl 

= (--)”m7T, 7Z = I, 2, . . . (33.19) 

The zero = o of the denominator is clearly exceptional 
since the numerator of (33.17) also vanishes for this value. 
The simplest way to see what happens in such a case, and 
to evaluate the corresponding term in the expansion formula, 
is to use the series for sinh ;2r and cosh z, namely 


cosh 2: == I H 7 1 


(33-20) 

(33-21) 


in (33.17), and cancel any common factor. Using (33.20) 
(33.17) becomes 

^_ ^pxlc)[x+{p^x^l6c^)+ . . .] 
p{pllc)[i + {pH^I6c^) 


£«[: -f (p^x^iec^) + • • -] 


(33-22) 


pl[i + (pVI6c^) +...]’ ^ 

and, in this form, the numerator of 'V does not vanish when 
^ = 0, and the denominator has only a single zero there. 




(33-23) 


Also 
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Using this result for ^ = o, and (33.19) for the zeros 
(33'I8), we get finally from (33.17) and the expansion 
formula (33.7) ^ 


V J- + 2 Conjugate^ 


Eaj 2E^( — )»» . tiTTX flrrct 
"7" i : > — ^ sin cos — r- - 


(33-24) 


The minor difficulty which was encountered above at 
/> == o arises in many problems, and it is always advisable 
to see precisely what happens at such a jpoint by substituting 
L (33*20), (33.:si), or the products (33.4), (33-5) for 

the hyperbolic functions, and cancelling any common 
factor. calculating the contribution from such a point, 
instead of using a term of the type r/(^r)^“’‘*/^^(^r)] which 
occurred in (33.7) and (2.8), one of the equivalent type 
(2.10), namely 


\ {p-ar)fm 
L g{p) 


^art 


(33-25) 


could have been used. In the problem above, this gives 


r Ej) smhpxlc l Ex 
Lp smhpl/c 

which is the first term of (33.24). 

Note also that if, after cancelling any common factors 
as in (33*22), the denominator of V has a multiple zero, 
a term of the type (2.20) must be used for this zero. 
Example 2, The submarine cable, e'=(^//c)^, where 

K — I /xvE. 

In this case (33.14) becomes 


(33-26) 


p sinh I y^(plK) ■ ■ 

^ in Example i, ^> = o is an exceptional value at which 
we denominator vanish. Using (33.20) 
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^ E«[l + (px^ISk) + ••■•] 

lp[i+{pPI6K) + . . 

-{-(pPISk) + . . .]}J g~ ^33-27) 

The other zeros of the denominator of (33.26) are 

p = — Kti^TT^jP, w = I, 2, 3 . . (33-28) 


Since 


sinh I 


\ 

Also [lli' 


± sinh ifiTT = o. 






2\k, 




(33.29) 


Thus we have from (33.26) and (33.7), using (33.27) and 

(33-29) 

V - + 5 ® T sinh«(j./K)* -| 

I ^ I liP I K)i cosh l{p I 




Ex 

T 




e-Kn^nnii^ sin 


(33-30) 


Example 3. The general case of (33.14) with q given by 
(33.15). Here = o is a zero of the denominator, and 


sii^h ql)~^ = [sinh ql]p^o == sinh Z(RG)^, ( 33 - 3 1 ) 

The other zeros of the denominator are found from the 
values of q which make sinh ql zero ; these are = o and 

^ = d: m == I, 2 . . . ( 33 * 32 ) 

The zero q 0 need not be considered, since a factor 
q has cancelled between numerator and denominator in 
(33.16). Using the value (33.15) of q, (33*32) gives 

(L^) + R)(C^ + G) + '^ = o, n = I, 2 . . ., (33-33) 
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This is a quadratic which gives two values of^ correspond- 
ing to each value of n ; it could also have been deduced 
from (33.16). Using the notation 

P - (R/2L) + (G/aC), a = (R/aL) ^ (G/aC), 

c = (LC)-*, . . ( 33 * 34 ) 

(33*33) becomes 

+ zpp + p2 — <7* + =0, n=i,2 (33*35) 

The roots of this are 

. . . (33 *,36) 

where v„ = [(«V®c®/Z^) — £r*] 4 , . . (33.37) 

and we suppose for simplicity that is teal for all values 
of n, that is ; if this is not the case (33.35) 

will have real negative roots for some low values oi n, and 
these must be treated separately. 

Hj, As usual we need 


[|(f siik 

.r«£+£lcoah,;l 


basing (33.31), we have finally by (33.7) 
V = E8inha!(RG)* 

smh/(RG)i 



{--Yine^vnt sin {mrxjl) 
Vn{— P + iv„) 



E sinh a;(RG)^ avrEc^ 
sinh /RG)* ~ir^ 


^ n(-)" 

;fev„(p^ + v„^)i 


sm COS (y^t — 


tanS„ = p/v„. . . (33.39) 


where 
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The examples above seem sufficient indication of the 
method of applying the expansion formula to problems on 
finite lines. From the results (33.13), (33.24), (33.30), 
(33.39) it appears that the expansion formula gives the 
solution as an infinite series, whose terms are products of 
trigonometric terms in x and exponential or trigonometric 
terms in t. These are quite different in form to those 
obtained in §§ 31, 3:^ in terms of multiply reflected waves. 
Each type of solution has its own advantages, so that it is 
hard to give any general indication of which method it is 
better to use — preferably both should be tried in any 
specific problem. However, the position as seen from the 
results of the examples worked out above is as follows : 

(i) For the lossless and distortionless lines with simple 
terminations, both methods are available, but the solution 
in terms of multiply reflected waves gives the actual wave 
form of the disturbance with a simple physical interpreta- 
ion, and thus is preferable. 

(ii) For more general lines the solution by the expansion 
formula is very little more complicated than that for the 
lossless line, while the solutions in terms, of multiply 
reflected waves involve much more complicated functions 
(not given here). 

.(iii) For lines with terminal impedances, the solution in 
terms of multiply reflected waves is useful in general only 
for the first few waves, that is, for not too large values of 
the time, while the solution obtained by the expansion 
formula (discussed in § 34) is available for all values of the 
time, though it is often difficult to determine the zeros of 
the denominator of its transform. 

34. The application of the expansion formula to the finite 

line with terminal impedances » 

Problems of this type are treated precisely as in §33, 
the only new feature being that the transforms have more 
complicated denominators whose zeros may be difficult to 
find. As a specific example we consider a lossless line of 
length Z, insulated at = /, and suppose that at ^ = o a 
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condenser of capacity Co, charged to voltage E, is discharged 
into the line at « = o. <= ’ e 

The subsidiary equation for the line is, writing c = (LC)-*, 

AiXT 


= o<x<l, 


with 


dx 


= o, X = 1 . 


(34-i) 

(34-») 


^ voltage and current at the point 

— 0 01 the line, the subsidiary equation for the condenser 




(34-3) 


Co is 

E 

that is, using (39.7), * 

^ “ f ’ • ^34-4) 

boundary condition to be satisfied at a; = o. 
ihe solution of (34.1) which satisfies (34.2) is 

? = A cosh [p(l - x)/c)]. 


Substituting this in (34.4) gives 


Thus 




EC„ 


coshp(Z — x)lc 


(34-5) 


^ ‘ ^h (pljc) + (pCJCc) cosh {pljcy 
The zeros of the denominator of (34.5) are ^ = o, and 

P==±tca.Jl, «=i, 3 (34.6) 

where the «„ are the positive roots of 

sin a + cos a = o, . . (34.7) 

A=Co//C . . . (34.8) 

IS the ratio of the terminal capacity Co to the capacity of 
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the whole line. Approximate values of these roots can 
be found from the intersections of the curves 

y = tan a\ 

and jy = — /eaj ’ 

as shown in Fig. 31, and these values can be improved by 
interpolation. 



Fig. 31 - 


To find V from (34.5) by the expansion formula (33.7), 
we need 


\^{ sinh ^ ^ cosh 

[_^dp C Oc C J Jp=icccn!i 

_ -I- cosh^ + 7^ sinh^l 

_^\cC cj c c Jp=.ica 

= {l/c)[k + I + k^oi^] cos a„, 
where k is given by (34-^)‘ 


• (34- 10) 
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Using {34.9) and (34.10), the expansion formula gives 


I k 


(*+i+^^a2)cosa, • ^ 


35. Further discussion of semi-infinite lines 

In § 30 the transforms of the solutions of problems on 
semi-infinite lines were found to contain the factor 

. . . ( 35-0 

where q = + ^)}> * ( 35 -^) 

and in §§31, 32 the travelling wave solutions for finite 
lines were found to involve similar factors. It is thus 
very important to know the functions which correspond 
to transforms of these types. In the case of the lossless 
and distortionless lines these caii be found by the use 01 
Theorem VI ; in all other cases, either the inversion 
theorem, or a table of transforms, must be used. Perhaps 
the most useful collection of transforms is given in Campbell 
and Foster^* Fourier Integrals for Practical Applications ’ , 
in the relevant part of which, Part 8, most of the results f 
can be read as Laplace transforms, p being our and ^ 
being written for our Many results for the general 
case (35.2) are given in their Nos. 860.5 onwards, and, 
for the submarine cable case, L = G = o, earlier. 

Here we consider only the submarine cable case, which 
corresponds also to linear flow of heat, and give an ex- 
tension of our table of transforms which will cover a 
number of problems of this type. The simplest and most 
important results follow from the known de^ite integral J 

f == . (SS'S) 

Jo ' 2,0 

* Bell System Technical Monograph, B-584. See also 
McLachlan and Humbert, Formulaire pour le calcul symbolique 
(Gauthier- Villars, 1941). 

t Namely, those stated for g > o. 
j Gibson, Treatise on the Calculust edn. 2, p. 470. 
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where a and h are positive. . Putting and sub- 

stituting u-=^ ^/tm this, gives 


Integrating (354) with respect to the parameter * 
from a to 00, gives 

foo foo 'n-i-f'® 1 

1 da — e~^^^da . (35.5) 

Jo Ja P^Ja 


Writing 


erf X 


2 

=— r e''‘^^du, 
TT-^Jo 


(35-6) 


2 

erfc :« = I — erf jc = — rl e~^^du, 


(3S-7) 


for the error function (numerical values of erf x are given 
in most books of tables) (35.5) becomes 




p-ap^ 


(35.8) 


In Table 2 the results (354) and (35.8) are given, together 
with some others which are useful in more complicated 
problems. Other transforms can be deduced from these 
by Theorems' I, III, IV, IX, also by differentiating or 
integrating with respect to a parameter, just as (35.8) was 
deduced from (354). For example (35.9) follows from 
(35.10) by differentiating with respect to the parameter a 
(or alternatively it may be derived independently from 
(35.3) by putting u = t'i'). Also (35.11) may be deduced 
from (35.9) by Theorem IV. The results (35.12) to 
(35.14) may be proved by a rather more elaborate discussion 
of the same type. 


^The processes of differentiating or integrating under the 
integral sign with respect to a parameter are most useful, since 
they often allow a chain of transforms (such as (3 5 - 9 ) ( 35 *ii)j 

or (35.12) to (35.14)) to be deduced from a single transform. 
It must always be understood that the use of these process^ 
needs justification, but this is usually readily supplied for integrals 
of the types considered here. 
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Finally, the Laplace transform of where v is not 
restricted to be a positive integer as in Table I, must be 
mentioned. This is 

. (3 s-i 6) 

Table z 


In this Table a and b are positive or zero 


y{p)‘ 

y[i). 



(35-9) 

p-ie-api 


(35-10) 

p-ig-Qpi 

erfc ( 

( 3 S-II) 

b ^ ' 

erfc -h bti^'j 

( 3 S-ia) 

g-ap^ 

p + bpi 

e«»+‘«erfc 

( 3 S-I 3 ) 

g-api 


( 35 -* 4 ) 

p(b p^) 

I 

pv 

fp-l 

wr • • 

(35-15) 


If v> — I, the integral in (35.16) exists and is denoted 
by P(v + i) ; this ^ves the result (35.15). Numerical 
values of the jT function are given in most books of tables - 
Its most important properties are 

r(v + 1) == vr(v),^ 

and r(n) — (« — i) !, if n is a positive integer, (35*17) 
so that (1.5) is a special case of (35.15), as it should be. 
Also 

Ai) = , and r(|) = . (35-1^) 

so that the special case a = o of (35.10) agrees with (35.15)* 
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As examples of the use of Table 2 we consider the follow- 
ing problems on the semi-infinite submarine cable. 

Example i. Constant voltage E applied at = o at 
t = o with zero initial conditions. Here, from (30.6) with 


L = G = 0, 

) 

1 

II 

(35-19) 

where 

k=i/RC. . 

(35-20) 

Then it follows from (35.11) that 




(35-21) 


Using (29.7), the current I is to be found frona 


R dx 


ECi 




(35-22) 


Thus, from (35.10), 

(l7Rf)t 

If we had only required the current 
becomes in this case 

T 


• (35-23) 

lo at * = 0, (35-22) 

• (35-24) 


and, from (35.15) and (35.18), 

'•-(S' • ■ 

Example 2. The problem of Example i, except that the 
voltage is applied to the line through an impedance 
As in § 30, the subsidiary equation for the line is 

= 0, x> o, . . (35-26) 

dx^ K 

where k is defined in (35-2°)., The solution of (35.26) 
which remains finite as oo is 

•? = . . (35-27) 
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Also, by (29.7), 

The constant A is found frona the boundary condition 

(32-2), 

Zal=^-'^, x = o. 


That is, 
Thus 


^Irk* + J p' 


■? = 


ERk* 




(35-29) 


+ R«*} 

If Zq is a resistance, R®, it follows from (35.14) ^l^^t 

( 35 - 3 °) 

where k = R/Rq. 

Again, if is a condenser of capacity Cq, (35*13) gives 

t 


' = E exp^ 


:} izVixt) + RCo/c*} 


R2CoV^RC(,/c 

(35-31) 

For other values of the expression i /(^op^ + R/c^) /n 
(35.29) must be expressed in partial fractions whose 
denominators are linear in^^. 


EXAMPLES ON CHAPTER IV 


I. Voltage E sin {oit + a) is applied at t = o through a ^re- 
sistance Ro to a semi-infinite lossless line « > o, with zero initial 
current and charge. Show that the voltage in the line is 


ERi 

Ro + Ri 


8m[a.(f-^)+a]H(r-^), 


where c = (LC)"*, R^ = V(L/C). 

2, Constant voltage E is applied at £ = o at the end x — o 
of a lossless line of length /, the end x " I being insulated. If 
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the initial current and charge are zero, show that the voltage 
at any point of the line is 




(zn -1- 2)1 — X 


7)- 


3* A lossless transmission line extends from a; ~ — / to jc = a, 
the end x — a being insulated, and the end x ~ — I being earthed. 
At i = o there is no current in the line, the portion o < x < a 
is charged to voltage E, and the portion — Z < k < o is uncharged. 
Show that the current at x ^ — I at time t is 


-® 7 (e)|.‘- ‘>-W' 


2n{l + a) + Z" 


H 




(zn -f i)(Z d- a) + a’ 


])■ 


where c — 

4. Constant voltage E is applied at £ = o through an inductive 
resistance Lq, Rq to the end jc = o o'f a lossless line of length Z, 
the end x — I being insulated, and there being zero initial current 
and charge. Show that the current at iw = o .is 


E 


Ro+Ri 


2ER1 


(Ro+Ri 




where c — (LC)“i, Ri = (L/C)i', and a = (Rq + Ri/Lp. 

5. A submarine cable (L = G — o, and k= i/RC in (29.10)) 
of length I has zero initial current and charge. The end jc = o 
is insulated, and a constant voltage E is applied at 5c = Z, show 
that the voltage at any point is 

E 4. y ±I_i)Lg-KC2n-l)w;4i« (an - i)^« 

TT ;J^(2M — i) zZ * 


6. A lossless line of length I has zero initial current and charge, 
and is insulated at ?c — Z. Voltage E sin a)£ is applied at a; == o 
for £ > o, show that the voltage at any point is 


E cos [co(Z — x)lc] sin oit 
cos 0)1 1 c 


00 

-}- SEcoZc ^ 

n—O 


sin [(z«+i)'n'c£/zZ] sin [(zn+ij^rJc/zZ] 
4co^Z® — (zn + ' 


9 
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for^any^w^ (LC)-i, and o) is not to be equal to (zn 4- i)‘jrclzl 

ar?A ^ submarine caMe of length / has the end x — I earthed, 
voltage E is applied at = o with 2ero initial gon- 
__ ■* using the method of § 31, that the current at 

X — O IS 
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